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A COMPUTER PROGRAM FOR HELICOPTER ROTOR 
KOISE USING LOWSON'S FORMULA IN THE 
TIME DOMAIN 


% 

C. Lucille Parks 
Langley Research Center 


SUMMARY 

A computer program (D3910) has been developed to calculate both the far- 
field and near-field acoustic pressure signature of a tilted rotor in hover or 
uniform forward speed. The analysis, carried out in the time domain, is based 
on Lowson's formulation of the acoustic field of a moving force. 

This report describes the digital computer program including the methods 
used in the calculations, a flow chart, program D3910 source listing, instruc- 
tions for the user, and two test cases with input and output listings ar.d 
output plots. 



INTRODUCTION 


When conditions of compactness of acoustic sources apply, the acoustic 
pressure signature of a rotor can he calculated by a fairly simple formula due 
to M. V. Lowson (ref. l). Although Low son 's formula has usually been used in 
the frequency domain, it can be used to compute a time history of the acoustic 
pressure signature due to steady and periodic forces on the rotor. 

This report describes a computer program developed to calculate both the 
far-field and near-field acoustic pressure signature of a tilted rotor in hover 
or uniform forward speed using Lowson *s formula in the time domain. Results 
obtained from this program agree with those available which use frequency 
analysis. Apart from the inclusion of forward speed, one advantage appearing 
from this program is its speed and simplicity. 

A description of the computer program for helicopter rotor noise is 
presented herein along with the methods used in the calculations, a flow chart, 
program listing, instructions for the user, and two test cases with input and 
output listings and output plots. 



SYMBOLS 


a vector whose components (a^, a^» a^) describe the location of the 
source point force on the rotor disk, meters (see fig. l) 

normalized amplitude of each harmonic obtained by the Fourier 
decomposition of the thrust 

phase angle of each harmonic obtained by the Fourier decomposition of 
the thrust, radians 

net force of the rotor blade on the air. Newtons 
velocity of the point force/speed of sound 

mach number of point force relative to the observer's direction 

✓ 

number of harmonics of the load 

2 

sound pressure, Hewtons/meter 

distance between the source point force and the observer, meters 
(see fig. 1) 

unit vector from the location of the source point force towards the 
observer 

effective radius of the rotor disk inscribed by the net force on the 
blades, meters (see fig. l) 

velocity of the rotor, revolutions/min. 


time the sound is observed, seconds 



T thrust, Hestons 

V. velocity of the vehicle, meters/sec. (see fig. l) 

n 

^ a vector whose components (x^» x^, x^) describe the position of the 

observer, meters (see fig. l) 

x q initial component of the observer's position, meters 

a ratio of the magnitude of the torque force to the thrust force 

# angle made by the normal to the rotor disk and the horizontal 

direction, radians (see fig. l) 

Q angular velocity of the rotor, radians/sec. (see fig. l) 

T time the sound is generated, seconds 

T initial generation time, seconds 

o 

At time increment 

6 blade angle, radians (see fig. l) 

6 initial blade angle, radians 

o 

A6 blade angle increment, radians 

Subscripts and superscripts: 

i denotes a computation for the ith source 

J denotes a computation for the Jth iteration 


k 




PROBLEM DESCRIPTION 

To apply Lowson's formula in the time domain for multiple sound source 
systems, that is, multiple blades on a given rotor, each source is followed 
separately and its acoustic pressure signature is calculated as a function of 
time. The overall pressure at any given time is the sum of the pressures which 
reach the observer at that time from the individual sources. Since there is a 
lag between the time that the sound is generated and the time that the sound 
arrives at the observer, care must be taken in executing this procedure. To 
get meaningful results, the arrival times of the acoustic pressure signals from 
all sources must overlap. 

In order to do this, the initial arrival time of the signal generated at 
the source nearest the observer must first be determined. Then each of the 
other sources is retraced back in time until its signal arrives within a 
specified range (determined by the number of harmonics and rotational velocity 
of the rotor) of the initial arrival time of the signal from the nearest source. 
This results in each source having its own initial generation time at which a 
signal is generated, and this signal arriving at the observer within a specified 
range of the initial arrival time of the signal from the nearest source. Thus, 
a time history of the acoustic pressure signature is computed for each source 
from its respective initial generation time. Then it is possible to sum the 
pressures from the individual sources to determine the overall pressure at any 
given arrival time. 

This application of Lowson's formula is described in reference 2. The 
computation of the components of Lowson's formula is described by the following 
equations : 


Blade Angle Increment 

if N = 0, A® 10 

if N > 0, A® (^) 


( 1 ) 


6 



Angular Velocity of the Rotor 

0 = 6 (sb) I 8 ™) 


( 2 ) 


Tima Increment 


At 


A© 

n 


( 5 ) 


After determining the initial generation time, t q , and the initial hlade angle, 
0 q , for each source, compute for each iteration, j, and each source, i: 


Blade Angle 


9? = 9 + 

1 o. 

i 


(j 


1) A© 


Generation Time 



l) At 


Observer Position (x^ component only) 



Source Location 



-R sin ©^ sin i 
R sin 0? cos § 
R cos 6? 

l 


GO 


( 5 ) 


(6) 


( 7 ) 
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Distance Between Source and Observer 

( 8 ) 


Observation Time 
r*J 

tj = (where 540 M/Sec. is the speed of sound) (9) 
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(14) 


( 15 ) 


(16) 
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Sound Pressure (Lowson's formula) 



Due to irregular arrival times of the signals, the stored values for each 
source must he interpolated at specified observer times and then summed to 
obtain the overall pressure signature. This is accomplished by using a 
second- order Lagrangian interpolation method. 

The calculations are continued for a specified number of iterations as 
explained in the Input Description. 
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PROGRAM DESCRIPTION 


The program D5910 is written in FORTRAN IV under the specifications of 
the Langley Research Center's SCOPE operating system for the Control Data 6000 
series. A FORTRAN variable list is given first, followed by a flow chart and 
a source listing. 

This program reads the input and determines the blade nearest the observer 
by calculating the minimum initial arrival time of the sound pressure signal. 
Each of the other blades is then rotated back in time until the initial arrival 
time of its signal is within a certain range of the minimum initial arrival 
time. Next, a time history of the sound pressure generated by eacn blade is 
calculated, starting with the respective initial arrival time of each blade. 

To determine the total sound pressure generated, the sound pressure for each 
blade is interpolated at specific time points using a second-order Lagrangian 
interpolation method. These interpolated values are then summed, giving the 
total sound pressure for the specific points. The program then writes the 
time history of the total sound pressure, and if the user so chooses, the 
program plots the time history of the total sound pressure . 

The interpolation subroutine, IUNI, and the plot subroutines, PSEUDO, 
CALPLT, AXES, LINE, and NFRAME, are called from the Langley subroutine library. 
Descriptions of these subroutines are included in the appendix, but the program 
listings for these subroutines are not included in this report. 



FORTRAH Variable List 

« 

A vector A, the location of the source point force on the rotor disc, 

in meters 

ALPHA a, the ratio of the magnitude of the torque force to the thrust force 

BETA B, a constant used to determine the number of sets of total sound 
pressure to he computed. (50 points per set) 

CGAMMA cos Y 

C0EFF1 C, the normalised amplitude of each harmonic obtained by the Fourier 
decomposition of the thrust 

C0EFF2 C* , the phase angle of each harmonic obtained by the Fourier decom- 
position of the thrust, in degrees 

CPHI cos f 

CSAV cos (iO + C'), i = 2, 3, . . N, N+l 

CTHETA cos 6 

DELTAU At, time increment , in seconds 

DELTET A8, blade angle increment, in degrees 

DIST r, the distance between the source point force and the observer, in 
meters 

DF 8F in Newtons/second 

3t 
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DM 

DOPFAC 

DIKES 

DTH 

ooa \ 

DUM2 » 
EMt 



PM 

GAMMA 

I 

IDHST 

IH» 

IPLOT 

IPOWER 


3M 

frr 

1 

1 - M • ? 

A0, blade angle increment , in radians 

dT 

dr 

Durany variables used in calling subroutine LISE 

M • r, Bach number of the source point force relative to the 
observer's direction 

M • M 

vector F, the net force of the rotor blade on the air, in Heatons 
P * M 

initial blade angle for each blade, in radians 
an index 

8 alphanumeric words of identifying information 
an error code word returned from subroutine IUNI 
an option to designate whether to plot 

used in determining the plotting scale factors for the X and I axes 


15 



”1 

J 

JPOVER 

K 


KC 

EJ 

KH 

EOUHT 

I | 

I 

HI 

SBLADE 

HL ) 
EP) 

HPTS 


OMEGAD 
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an initializing indicator used in calling subroutine IUNI 

used in determining the plotting scale factors for the X and Y axes 

used to count the number of points computed per blade 

used in determining the minimum values for the X and Y axes 

the number of sets of total sound pressure to be computed (50 points 
per set) 

used to count the umber of points computed per blade 

an array containing the number of points computed for each blade 
indices 

the number of harmonics of the load 

H ♦ 1 

the number of blades on the rotor 
indices 

the number of points to be plotted 

used to save the last 5 points for the next interpolation set 
Q, angular velocity of the rotor, in degrees/second 



2 

sound pressure, in Hevtons/meter 


OMEGAR 

P 

PHI 

PHIR 

PMAX 

BOH 

POFT 

PTAB 

RADIUS 

RHAT 

RHDF 

RHDM 

RHF 

* 

ROC 


< , the angle made by the normal to the rotor dish and the horizontal 
direction, in degrees 

i , in radians 

the maximum interpolated value of total sound pressure, 

2 

in Hevtons/meter 

the minimum interpolated value of total sound pressure, 

2 

in Hevtons/meter 

the interpolated value of sound pressure returned from subroutine 

2 

IUHI, in Hevtons/meter 


a table used to store values of sound pressure for interpolation 
and for plotting 


R, the effective radius of the rotor disk inscribed by the net 
force on the blades, in meters 


r, a unit vector from the location of the source point force towards 
the observer 



£ • P 

BO 

340 . 


in radians 



HOT 2 

HOSC — SUSP ^ radians /se: 

RFN velocity of the rotor, in rev/min 

SAVE iB + C{, i = 2 , 3, . • H, H+l 

SECTOR — — in radians 

HBLADE 

SGAMCA sin Y 

SPHI sin f 

SSAV sin (i6 + C'), i = 2, 3, . . N, N+l 

STHETA sin 6 

SIML cos # + a cos 6 sin f 

SUM2 -sin i + a cos 6 cos 9 
SUM3 a ft T sin 0 

T t, time the sound is observed, in seconds 

TAU T, time the sound is generated, in seconds 

TEST initial T for each blade, in seconds 
TH t, thrust, in Hevtons 

THETA 6, blade angle, in radians 

TIME time at which the sound pressure is interpolated, in seconds 
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TMAJ 

TMAX 

mra 

BOR 

TO 

TOTAL 

TTAB 

VH 

VECM 

VHC 

X 

xo 

XSF 

YSP 


distance between major tic marks on the plotting axes* in inches 

used as maximum time for computing sound pressure and as maxim® 
interpolation time* in seconds 

initial arrival time of the sound produced by the blade nearest the 
observer* in seconds 

number of divisions per inch for minor tic marks on the plotting axes 

initial interpolation time point for each interpolation set* in 
seconds 

an array of the summed interpolated values for all blades at specific 
time points 

a table used to store values of arrival time for interpolation and 
for plotting* in seconds 

V. , velocity of vehicle, in meters/second 
n 

M, velocity of the point force/speed of sound 



340 . 


x 1 * Xg, x^, the position of the observer at any given time T, in meters 
x q , the initial component, in meters 
scale factor for the X axis 
scale factor for the f axis 
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PROGRAM ROTOR I INPUT* OUTPUT* TAPE5« INPUT, TAPE6-0UTPUT , TAPEIOI 
C 

DIMENSION AI3I * COEfFl(SO)* C0EFF2I50I, OfU). 0MI3I* FI3I* 

IGAMMA (41* IOENTIBl* FOUNT (4 I, PI200*4>, PTAB(202*tl* RMAT<3>* 

2T( 200*41* TESTI4I * TOTAL I 2001* TTAB<202>* »ECN(i). X«3» 

C 

C 

NAMELIST /OAT A/ ALPHA* BETA* COEFF1. C0EFF2* IPLOT* N* N0LAOE* 

1PHI* RADIUS* RPM* VH* X 

C 

C **• INPUT **♦ 

C IC0EFF2 MUST BE GIVEN IN DEGREES. THESE ARE THE PHASE 

C ANGLES OF EACH HARMONIC STARTING FROM FUNDAMENTAL. 

C SET THE PHASE OF FUNDAMENTAL TO 0.0) 

C 

CALL PSEUOO 

500 REAO 15,501) IDENT 

501 FORMAT (8A10I 

IF IE0F.5I 502*503 

502 IF IIPLOT.EO.il CALL CALPtTI0.*0.*999» PLOT 

STOP 

503 REAO (5.0ATAI 
WRITE 16 .OATAI 
PRINT 504* IOENT 

504 F0RMAT(1H1*2X*8A10/// I 
C 

C •«* INITIAL I /AT ION 

C 

XO ■ XII) 

EN * N 

IFIN.EQ.O! EN*1» 

DELTET * 10. / EN 
N1 * N»1 

0 ME GAO * RPM * 6. 

DELTAU - OELTET / OMEGAO 
OTETR = OELTET * .0174532025 
OMEGAO » OMEGAO • .0174532925 
PHIR » PHI * .0174532925 
CPHI » COSI PHIRI 
SPHI « SINIPHIRI 
DC 10 I-1.N1 

10 C0EFF2I T I - C0EFF21 1 1 * .0174532925 
VHC » VH / 343. 

ROC * I RADIUS • OMEGARI / 340. 

RCSC » ROC * OMEGAR 

SECTOR * 6.283105306 / NBLAOE 

IF (IPLOT.Ed.OI GO TO 40 


PMIN * 1.E1G PLOT 

PMAX * l.E-10 PLOT 

NPTS * 0 p LOT 

REMINO 10 PLOT 

C 


C **«• DETERMINE INITIAL TIME POINT FOR INTERPOLATION *«* 

C 

40 DO 50 L*l»NBLAUt 

GAMMAIL I - IL-ll • SECTOR 
C GAMMA - CJS (GAMMA! LI I 
SGAMMA - SINIGAMMAUI) 

Alii « -RADIUS * SGAMMA • SPHl 
A12> » RADIUS * SGAMMA * CPHI 
AI3I * -RAOIUS ♦ CGAMMA 

DIST » ( IX0-AI1I |**2 ♦ 4X12 I-AI2) 1**2 ♦ U13»-AI3II*«2»»».5 
TESTILI • 01 ST/ 340. 

50 CONTINUE 

TNIN * 1.E10 
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00 51 L*1 .NBLAOE 

51 IF itestili.lt.tnini tnin-testili 

*** ADJUST TAU FOR EACH BLAOE **• 

00 75 L"*l» NBLAOE 
1*1 

tau * 0. 

GANMAIL) * IL-ll • SECTOR 
60 IFIITEST(L)-TMIN).LE.OELTAU) GO TO 74 
IFII.LE.1000) GO TO 70 
PRINT 505.1 

505 FORMAT! 1H1.*TIM£ HAS HOT REACHED TMIK FOR BLADE* 12) 

STOP 

70 TAU * TAU * OELTAU 
Kill • M - « • TAU 
GANMAIL) * GANMAIL) - OTETR 
C GAMMA * COSIGAMMA(O) 

SGAMMA * SINIGAHMAILII 
All I « -RADIUS * SGAMMA • SPH1 
AI2) • RADIUS • SGAMMA • CPHI 
AI3I * -RADIUS • CGAHMA 

DIST * (U(l>-Atl))**2 ♦ IX(2I-A( 2) )**2 * IXt 3)-AI3)l**2» **.5 
TESTIlI - TAU ♦ 01ST/340. 

1*1*1 
GO TO 60 

74 TESTIl* * TAU 

75 CONTINUE 

*** K IS THE NUMBER OF SETS OF INTERPOLATIONS *** 

PRINT 507 

507 FORMAT!// /7X.*T!ME*15x.*tNTERP0LATE0 PRESSURE SUN*) 

AC * o’ 

IF IVH.EQ.J.O) GO TO 100 
A = 1FIXIBETA*X0 / IVH*0£LTAU*50>> ♦ 1 
GO TO 101 

100 A a IFIX 1360 / (NBLAOE • DELTET • SOI! * 1 

101 00 1002 M * ItA 
THAT = THIN ♦ I49+IH-1 )*50)*DELTAU 

•*• NBLAOE IS THE NUMBER OF BLADES *•* 

00 104 L*l. NBLAOE 
J 3 1 

*** BLADE ANGLE. RELATIVE TINE. A NO Xtl) COMPONENT OF OBSERVER *** 

. 

105 THETA - GANMAIL) * IJ-1) * OTETR 
CTHETA * COStTHETAI 
STHETA * SINITHETA) 

TAU * TtSUL) ♦ IJ-1) • OELTAU 
XII) = XU - Vrt * TAU 

•*• A -VECTOR, OISTANCE. ANO ABSOLUTE TIME 

All) * —RADIUS • STHETA • SPH1 
A(2) * RADIUS • STHETA • «J>Ht 
Al 3) * -RADIUS • CTHETA 

01ST » 1 1 XI 1 )-A( 1) |**2 ♦ IX(2)-AI2))**2 * IXC3)-AI3>)**2) **.5 
KJ * KC ♦ J 

TIKJ.L) • TAU *• 01 ST / 340. 

C 

c •*• N-VECTOfl. RH AT- VECTOR » ANO THEIR OOT PRODUCT **• 

C 
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veCMUJ = VH/340. - ROC * CTHETA * SPHI 
VECMI2) * RCC * CTHETA * CP HI 
VECNI3I » ROC • STHETA 

ErtSO = VECNm*VECMfl) ♦ VICHI2i»«ECHl2l ♦ VECMf 3)*VECMI3) 

c 

RHAT 111 » 1X11) - Alin / dist 

RHATI2I * 1X121 - AI2II / DIST 

RHATI3I = ( XI 3) - AC 31 1 / DIST 

C 

ERR - VECM 1 ) *RHAT I II ♦ VECMt 2I*RHATI2> » VECM(3I*RMAT<3) 

IF (EMR.LT. l.OI GO TO 103 
PRINT 102 

102 FORMAT (* MACH NUMBER (S NOT LESS THAN ONE. PROGRAM ABORTED.*) 

STOP 

*** PARTIAL DERIVATIVE OF M WITH RESPECT TO TAJ **• 

103 OMIl) = ROSC * STHETA * SPHI 
0MI2) * —R3SC * STHETA * CPHI 
DM3 ) * ROSC • CTHETA 

THRUST AND ITS DERIVATIVE 

TH = CQEFF1I1I 
OTH * 0. 

IF (N.EU.OI GO TO 110 
00 104 1=2, N1 

SAVE = I • THETA ♦ C0EFF2C I ) 

SSAV = SINI SAVE I 
CSAV * COSC SAVE) 

TH = TH ♦ C0EFF1II) * CSAV 

104 DTM a DTH - I • C0EFF1I1I • SSAV • OMEGAR 

*»* F-VECTOR *** 

110 SCJMl = CPHI #• alpha*ctheta*sphi 
SLM 2 * -SPHI * ALPHA*CTHETA*CPHI 
Fill a -TH • SUM1 
FC2) * TH • S0M2 
F( 31 a ALPHA * TH * STHETA 

FM = F(1)*VECM(1) » FI2 )*VECMI2 ) ♦ F(3)*VECM(3> 

*** PARTIAL 0EP1VATIVE OF F WITH RESPECT TO TAU •*• 

SOM 3 = ALPHA • OMEGAR • TH • STHETA 
OF 111 » -DTH • SUM1 ♦ SUM3 * SPHI 
L'F 1 2 1 a OTH • SUH2 - SOM 3 • CPHI 

OF 1 31 = ALPHA * OMEGAR * TH * CTHETA ♦ ALPHA * DTM • STHETA 
SCUND PRESSURE **• 

RhOF a RHATIll* OF(l) ♦ RHATI2I*0F<2) ♦ RHATI3)*OF(3) 

OOPFAC a 1. / (l.-EMKI 

RHF = ( RHAT I 1 l*FI 1 ) *■ RHA T|2)*FI2» RhAT|3)*F|3>| • OOPFAC 

RrtOM a RHAT(1)*DM(1| ♦ RHAT(2)*CM(2) ♦ RHAT 1 3) *0M( 3 1 
PCKJ.LI = .074577471© * I OOPFAC** 2/DI ST ) • (RFJF/3A0. *■ RHF*RHOM 
1/340. * (ll.-EMSU) * RHF - FM) / DIST) 

IF (ITMAX - TIJ.LII .LE. DEL TAU) GC TO 108 
IFIKJ.GE .200 ) GU TO 106 
J»J*1 
GC TO 105 

106 PRINT 107 ,L«M 

107 FORMAT (LH1,*BLA )E *I1,« FAILED TO REACH TMAX iN 200 ITERATIONS FT 
IK INTERPOLATION SET*, 131 

STOP 
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108 GAMMA (LI « THETA ♦ OTETR 
TESTILI » TAU OELTAU 
KOUNTILI * KJ 

109 CONTINUE 

**• INTERPOLATION *** 

00 234 NP=1 .50 
204 TOTAL INPI * 0. 

TO » THIN ♦ (M~1I*S0*0ELTAU 

00 220 NL*1 .NBLAOE 

TIME * TO 

KN * KOUNTINLI 

DC 206 NP«l,KN 

TTAtHNP) = T(NP.NLI 

206 PTAdINP.il * P(NP,NLI 
IPT a -1 

OC 220 NP-1.50 

207 CALL IUNI (200.K0UNT (NL I «TTA8. l.PTAB. 2. TIME.POFT .IPT.1ERRI 

*•* CHECK ERROR COOE •** 

IFI IEAR.E0.0I GO TO 213 
Ir(IERR.LT.O) GC TO 207 
PRINT 208* I ERR 

200 FORMAT (1H3,I5.*Th element out of oroer. no (nterpolation.«i 
STOP 

209 IFRR a IABSIIEHRI 

GO TO (210. 211. 212.2131. tERR 

210 PRINT 214 

214 FCRMAT (1H).*IOK0Ek ■ 0*1 
STOP 

211 PRINT 215 

215 FORMAT PlHO.*ONE ELEMENT IN X ARRAY* I 
STOP 

212 PRINT 216 

21b FORMAT (1H0. * INSUFF IC IENT NUMBER OF POINTS. NO INTERPOLATION.*! 

STOP 

213 TOTAL (NP) * TOTALINPI ♦ POFT 
IFINL.LT.N3LA0E) GO TO 217 
PRINT 5Jd. TIME. TOTALINPI 

SOB FORM ATI 2X .E 15.7 . 14X.E15 .7 1 
IF IIPLOT.EO.OI GO TO 217 


MRITE (101 TIME, TOTAL INPI PLOT 

NPTS » NPTS ♦ 1 PLOT 

tFITOTALINPl.LT. PMINI PMINaTOTAL (NP I PLOT 

IF! TOTAL (NPI.GT.PMAXI PMAX* TOTAL INPI PLOT 

217 TIME » TIME ♦ OELTAU 
220 CCNTINUE 


*** SAVE LAST 5 POINTS FOR NEXT SET OF TABLES **• 

230 OC 25 J NL*l*NdLAOE 
NS a KOUNTINLI -5 
OC 250 NPal.S 
TINP.NLI a T(NP»NS.NLI 
250 P (NP.NL I » PINPFNS.NLI 
KC * 5 

1002 CCNTINUE 

IF (IPLOT.EO.OI GO TO 500 


C PLOT 

C ••* PLOT TIME VS. SOUNO PRESSURE *** PLOT 

C PLOT 

TMaX a TIME - OELTAU PLOT 
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END FILE 10 
REMIND 10 

*•* DETERMINE scale factors *** 

XSF « (TMAX-TMIN 1 / 22 . 

IPOMER » 0 

1005 IFIXSF.GT.. 9999999901 GO TO 1006 
XSF • XSF*lO. 

IPOMtR a IPOMER * 1 
GC TO 1005 

1006 IX > IFIXIXSFf.il 

XSF a FLOAT ( IX)/10.**IP0MER 
JPOmER a o 

1065 IF ( JPOMER. GE. IPOMER) GO TO 1066 
THIN a TM1N*10. 

JPOmER > JPOMER ♦ 1 
GO TO 1065 

1066 IX a IF IXt TMINI 

TPIN - FLOAT UX)/10**JP0mER 

VSF a (PMAX-PMINI/9. 

IPOMER a o 

1007 IFIVSF.GT.. 999999991 GO TO 1000 
VSF a VSF * 1J 

IPOMcR a IPOMER ♦ 1 
GC TO 1007 

1003 IV a IF ixl VSF*. 51 

VSF « FLOAT I IVI/10.** IPOMER 
JPOMER a 0 

1065 IF(JPGmER.GE.)POmEK) GO TO 1006 
PMIN a PHI N* 10* 

JPOmER a JPOmER ♦ 1 
GO TO 1005 

1086 IV a IFIX(PMIN~.S) 

PMIN a FLOAT I IV l/10«*JP0MER 

*** ORAM AXES 

TMAJ a 1. 

TMIR a U. 

CALL CAL PL T I1..1..-3I 

CALL AXES (0.t0..0.0.25..TMlN.XSF.TMAJ»TMIR,AHTIME* *2.«A) 
CALL AXES <O*»O.»9O.*lO.»PMlN»VSF ( TMAJtTMIR*8HPRESSURE'.2f0) 

*** REA C ANl) PLOT POINTS •»* 

J * 0 

00 2005 lal.NPTS 
J « J ♦ 1 

READ (10) TTAB(J). PTA6IJ) 

IF( I.GE.NPTS) GO TO 2000 
IF(J.LT.200) GO TO 2005 
2000 TTA6( J»l) a THIN 
PTAdlJt-l) a PM|N 
TT AB( J »2I * XSF 
PTAb ( J*2 I * VSF 

CALL LINE (TTAB, PTAtt * J .1 ,0 .DO**, 00M2 I 
TTAB(l) a TTAtfIJI 
PTAB(l) a PTABIJ) 

J * 1 

2005 CONTINUE 

CALL NFRAME 

GC TO 5 JO 
END 
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PLOT 

PLOT 

PLOT 

PLOT 

PLOT 

PLOT 

PLOT 

PLOT 

PLOT 

PLOT 

PLOT 

PLOT 

pljt 

PLOT 

plot 
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PLOT 

plot 

plot 
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PLOT 
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USAGE 

Program Information 

The program D39IO is run on the Control Data 6000 series under the SCOPE 
operating system. If the plot output is chosen, a 10” x 25" plot of the total 
sound pressure time history will be executed on the Varian plotting device. 
Minimum machine requirements are 550000 locations of core storage. The 
time required to compile the program is approximately 2.0 seconds. If the plot 
option is not activated, approximately .02 seconds sure required to compute one 
time point of total sound pressure, and one 0/S call is required for approxi- 
mately 6 points. If the plot ou.prt is chosen, approximately 2 0/S calls sure 
required for each point, and the time required for computing per point increases 
slightly. 
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Deck Configuration 

The following sketch shows the deck configuration needed for execution: 
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The following sketches show the control card configuration needed for execution both with and 
without plotting: 


with plotting: 


without plotting: 



'PLOT VARIAN 



Igo. 


LJ 

1 (RUN card 

-j 


PLOT. VAR IAN 


RFL( 52000) 


LGO. 


f USER card 


JOB card 


^ RUN card 

'norfl. ~~ 

fuSER card I 


(JOB card 
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Input Description 

The first card of input for each case is 80 columns of alphanumeric data 
which is printed out as an identification heading. 

The succeeding cards of input use the FORTRAN NAMELIST capability with 
DATA as the NAMELIST name. The maximum allowable dimension appears following 
the variable name, unless it is a single value variable. Minimum and maximum 
allowable values and required units of the variable name are also given, if 
appli cable. 

ALPHA the ratio of the magnitude of the torque force to the thrust force 

RADIUS the effective radius of the rotor disk inscribed by the net force on 

the blades, in meters 

NBLADE number of blades on the rotor, 1 £ NBLADE U 

PHI angle made by the normal to the rotor disk and the horizontal direc- 

tion, in degrees 

RPM velocity of the rotor, in revolutions/minute 

VH velocity of the vehicle, in meters /second 

X(3) initial position of the observer, in meters 

N number of harmonics of the load, 0 <_ N £ U9 

BETA a constant used to determine the number of sets of total sound 

pressure to be computed (50 points per set.); used only if VH = 0 

NOTE: K = IFIX (BETA * X0/(VH*DELTAU*50) ) + 1 

where K = number of sets 

XO * initial X(l) component 

DELTAU = 10 

N*RPM*60 
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C0EFF1(50) normalised amplitude of each harmonic obtained by the Fourier 
decomposition of the thrust; N+l amplitudes must be supplied 

COEFF2(50) phase angle of each harmonic ob' '.ined by the Fourier decomposition 
of the thrust; in degrees; N+l phase angles must be supplied (set 
the phase angle of the fundamental to 0.0) 

IPLOT plot option 

IPLOr = 0; no plot 

IPLOT 4 0; plot sound pressure vs. time 
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Output Description 

The output of program D3910 consists of printing and optional plotting. 

For each case, the BAMELIST input is printed first. Then the alphanumeric 
identification input is printed, followed by a time history of total sound 
pressure. The time points that are printed are the times at which interpolation 
occurred. The total sound pressure that is printed is the sum of the inter- 
polated values of sound pressure of all blades. 

If IPLOT = 0, no plotting will occur. If IPLOT 4 0, a 10" x 25" Varian 

plot representing the time history of total sound pressure will be part of the 
output. 
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Input for Sample Case 1 


SAMPLE CASE 1 — VH * IS* 

tUATA ALPHA * .114, RAUlUS * 4*3* ft BLADE * 4* PHI * i6.t RPM * 355*. 
VH » IS* * X * 30* • 5 0* , 10* • ft * 15* BETA * *25* IPLOT » 1* 

COcFFl * 1.. *4o236, .0537C* .05012* .03273* .01905, .00031* .01585* 
.02461. .10708, .00794* .00631, *01776* *00531. *00501. *00467, 3**0*. 
C06FF2 * 4*. 40., -25* * 20*. 15* • -30* , 44*0.* 


Output for Sample Case 1 


SCAT a 






ALPHA * 

0 . 1 I 4 F* 00 , 





*CTA 

O.25E*O0. 





COEFFl - 

0.1 c *0l, 0. 

46238 c *00» 

0.53? e — Cl * 0.5C12F— 01* 

0.32I3E— 01* 


0.1405E-01, 

0.841E-02. 

0. IS85E-01. 0.2661F— 01 

* 0.1C8E-C2 

* 


0.79AE-02, 

0.63IE-02. 

0.i77eF-01, 0.5315-02, 

0.5GIE-C2* 



0.46IE-C2. 

0*0* 0 .0 9 

0.0* O.C, 0.0, 0.0, 

0.0* 0.0* 

0.0* 


0*0« 0#0» 

0*0t 0.0* 

0.0* 0.0, 0.0, 0.0, 

0.0* c.c * 

0.0* 


0.0* 0.0* 

0.0* O.Ot 

4.0, O.C* O.C, C.O, 

0 .0 * O.C* 

0.0* 


0*0* 0*0* 

0.0* 




CnPFF 2 s 

0.0. 0.4r*02, -0.25 c *02» 0.25*02, 0.15 C *02, - 

C.3E+02* 0. 

0* 

o.r, o.o. 

0.0. o.o. 

0*0* 0.0* 0.0* 0.0# 

0 .0 * O.C* 

0.0* 


0.0* 0 • 0 * 

0.0, 0.0, 

0.0* 0.0* O.C* 0.0* 

C.O* O.C* 

0.0* 


C.O* 0 - 0 * 

0.0, 0.0. 

0.0* 0.0* 0.0* 0.0* 

0.0* O.C* 

0.0* 


0 *0 * 0.0 * 

0.0, 0.0, 

0.0* 0.0* 0.0* 0.0* 

C.O* 



IPLOT * 1, 

ft * IS. 

NHLAOE * 4* 

PHI = 0.26 c *02* 

RADIUS = 0.43E*Ol, 

RPM = 0.3SSF+03, 

VH - 0.155*02, 

% - 0.3F *02 * 0.5» : »02, 0.1F*02, 

tFftD 
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SAMPLE CASE I - VH * 15. 


TIME 

INTERPOLATED PRESSURE SUP 

I.6754662E— 01 

-1.10482036-04 

1.67859616-01 

-1.00562546-04 

l .681 7260f— 01 

-9. 3 596051 6-05 

1.68685596—01 

-8.97551 75*-05 

1.68798586-01 

-8.96118166-05 

1.69IH57 C -01 

-9.35478786-05 

1. 6942.56 E-01 

-1.01817266-04 

l .69 73755 E-01 


1 . T005054E-01 

-1.31 709056-04 

!. 70363536-01 

-I. 52996486-04 

1.706T651 C — 01 

-1.78007726-04 

1 • 7098950E -01 

-2.00025526-04 

1.71 302696-01 

-2.3622509P-04 

1. 71615686-01 

-2.676895C6-04 

1.71928676-01 

-2.9 9150386-0* 

1. 7226166E-01 

-3.29440676-04 

1 • 72556656-01 

- 3. 57359426—04 

1.7286766^-01 

-3.8148945E-04 

1.73180636-01 

-4. 007002 36- 04 

1.73693626-01 

-4.1395143P-04 

1.73806606-01 

-4 .20201566—04 

1.76119396-01 

-4.18802566-04 

1 • 7663238E-01 

-4.09316726-04 

1. 7676537F-01 

-3.91 58525E-04 

1 . 1505836 --01 

-3.65171296-04 

1.7537135E-01 

-3.3229565E-04 

1.7568636=-01 

-2.9203063 c — 04 

1.75997336-01 

-2. 4 5986 096-04 

1.76310326-01 

- 1.953869*6-04 

1. 16623316-01 

-1. 41891996-04 

1.76936306-01 

- 8. 6942 9 6 5E— 05 

1.1726928 c -01 

-3.23940846-05 

1.77562276-01 

2. 01542556-05 

1.77875266-01 

6.9214152E-05 

1.78188256-01 

l. I321501F-04 

1.78501266-01 

1.51051026-04 

1.78816236-01 

1.81826996-04 

1.79127226-01 

2. 0488657 6—04 

1. 7966 02 1 6-01 

2.1 99909 75— C4 

1.7975320E-01 

2 .2 72 5 5 7 6 6-04 

1.80066196-01 

2.2709674F-04 

1.80379176-01 

2 ♦ 201 8936 E- 04 

1.80692 I 66-01 

2. 076 361 OF— 04 

1.81005156-01 

1.90497916-04 

1 . 81 31 816 c -0l 

1.70275636-04 

1.8l63113 c — 01 

1.40292876-04 

1.81066126-01 

1.25878596-04 

1.8225711-— 01 

1.044941 7E- 04 

1.82570106-01 

0.52 108686—05 

1.82883096-01 

6.66723006-05 

1.83196086-01 

5. 66757596-05 

1.83509)76-01 

4.06245366-05 

1.83822056-01 

4.50626986-05 
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1.841 35046-01 
1. 8444803 c -0 l 
1.84761026-01 
1.85074016-01 
1.85387006-01 
1 .85699996—01 
1.8*01298^-01 
I • 8*325976— 01 
1 • 866389* e -0l 
1.86951946-01 
1.87264936-01 
1.87577926-01 
1.87890916-01 
1.88203906-01 
1.8851689F— 01 
1 • 8802580*— 01 
1.89142876-01 
1.89455866-01 
1.89768856-01 
1.90081846-01 
1.50394 82*-01 
1.5C707816-01 
1.9IC2C80 c -0l 
1.91?? 3796-01 
1.91646786-01 
1 .91959776-01 
1 .92? 7? 766-01 
1.92585 75 c -01 
1.92898746-01 
1.93211736-01 
1.93524716-01 
1.03837706-01 
1.94150696-01 
1.94463686-01 
1 .94 7 766 76-31 
1 . 95^89666-01 
1.95402656 -01 
1.95715646-01 
1.96028636-01 
1.96341626-01 
1.96'- 54(16-01 
1 .9696 7S9--01 
1.97280586-01 
1.97593576-01 
1.97906566-01 
1.98219556-01 
1.98632546-31 
1.9884553^-01 
1.9915852F-01 
1.99471816-01 
1.99784506-01 
2.000 17886—01 
2.00410476-01 
2.00723466-01 
2.01036456-01 
2.01349448-01 
2. 31662436-01 
2.01975426-01 
2.02288416-01 
2.02601606—01 
2.02914396-01 
2.03227386-01 


4.59035966-05 
5.07779926-05 
5.91379856-05 
7.02897326-05 
8. 3>9071 3F— 05 
9.73137276-05 
1. 1 133336E— 04 
1.24519386-04 
1.36111*0*— 04 
1.45390726-04 
1.5191209^-04 
1.55346856-04 
1.55575786-04 
1.52638416-04 
1.46773646-04 
1.382876*6-04 
1.27692926-04 
I. 15466356-04 
1.02127826-04 
8.828465*6-05 
7.43923866-05 
6.08976166-05 
4.81618046-05 
3.64108646-05 
2.57996766-05 
1.63544636-05 
8.03207286-06 
6.43894766-07 
-5.87569636-06 
-1.20719006-05 
-l. 79571806-05 
-2. 38 779716-05 
-3.00657286-08 
-3.67708436-05 
-4.41302266-05 
-5.22358876-05 
-6.11433986-05 
- 7. C 3 108806— 05 
-8.11658956-05 
-9.20614516-05 
-l.033264»6-04 
-1.14768146-04 
-1.26153096-04 
-1.37286226-04 
-1.47953236-04 
-1.58018216-04 
-1.6/308906-04 
-1.75743236-04 
-1.83280526-04 
-1.89903986-04 
-1.95654456-04 
-2.00607706-04 
-2.04864116-04 
-2.00543066-04 
-2. 1 1 7693C6— 04 
-2.1466664C-04 
-2.17360846-04 
-2.19935706-04 
-2.22442016-04 
-2.24899366-04 
-2.2 7282996-04 
-2.29532686-04 



2. 03540366-01 
2.03853356—01 
2.0416614F-01 
2.044 79336-01 
2.0479232F-01 
2.05105316-01 
2. 05410306-01 
2.05731296-01 
2. 060442 8F-01 
2.06357276-01 
2.066 70256-01 
2. 06983246-01 
2.07296236-01 
2. 07609226-01 
2. 079222 1F-01 
2. 08235206-01 
2. 08 548 196-01 
2*08861 1 8*:— O I 
2.09174176-01 
2.09407166-01 
2 .C980C156— 01 
2 . 101 131 36-01 
2.13426126-01 
2 ♦ 107391 l F— 01 
2 .11052106-0 I 
2.11 36509F— 01 
2.116 78086-01 
2.11991076-01 
2.12304066-01 
2. 126 1 7056-01 
2.12930046-0! 
2.13243026-01 
2. 13556016-01 
2 . 1 3 86900 E- 01 
2 • 141 81 576-01 
2.1649498F-01 
2.14807976-01 
2.15! 20966-01 
2. 1 54 3 395 c — 01 
2.15746946-01 
2. 160599 36-01 
2. 163 12926—01 
2.16 6 85 906 -91 
2.l699889 r -01 
2.17311886-01 
2. 176248 7F-01 
2.! 793786 F-01 
2. 182 5085 F-Oi 
2.18563846-01 
2.18076836-01 
2.1918982^-01 
2.19502816-01 
2.I981579F-01 
2.20128786-01 
2.20441 776-01 
2.20754 766-01 
2.21067756-01 
2.21380746-01 
2 • 21 693 73E-01 
2.22006726-01 
2.22319716-01 
2.22632 706-01 


-2.31519236-04 
-2. 33115786-04 
-2.34124626-04 
-2.34355366-04 
- 2. 336222 6F— 04 
-2.31671716-04 
-2.28375 386—04 
-2.23587406-04 
-2.17211186-04 
-2.09232616-04 
-1.99704396-04 
-1.88770176-04 
-1.766051 7F-04 
-1.63749076-04 
- 1. 60421276— 04 

- 1.37182266-04 
-1.24589646-04 

- 1.13213205— C4 
-1.03758216-04 
-9.68188006—05 
-9.29621 59F-C5 
- 9*26897556- C 5 
-9.64 748326-05 
-1.04573 746-04 
-1.17118576-04 
-1.34074406-04 

- 1.55192896-04 
-1.80081246-04 
-2.0 8016296-04 
-2.3819^006-04 
-2.69658196-04 
-3.01292216-04 
-3.31740676-04 

- 3. 5981615^—04 
-3.34255356-04 
-4.03721066-04 
-4.1 7195166-04 
-4.23650756-04 
-4.22 39 7866-04 
-4.12971386-04 

- 3. 95 19 70 8F— 04 
-3.69226996-04 
-3.35472576-04 
-2.94806346-04 
-2.48253006-04 

- 1.97032636-04 
-i. 42853916-04 
-8.71787726-05 
-3.18908916-05 

2.13807166-05 
7.11122506-05 
1.15700386-04 
1.54C09736-04 
1.05127466-04 
2.08381066-04 
2.23528016-04 
2 .306077 8 6-04 
2 .30289206-04 
2.23028296-04 
2.10033346-04 
1.92399356-04 
1.71663346-04 
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2.22945698-01 
2. 252566 78-01 
2*23 571666-01 
2* 23884656—01 
2.241 9764 6-01 
2.2451063 c -01 
2 • 2482 362 E*0 1 
2.25136616-01 
2. 2544 96 0€ *01 
2 .25762596—01 
2. 260755 8E-01 
2.26388568-01 
2.26 701 55 c — 01 
2.27014546-01 
2.27327536-01 
2.27640528*01 
2.27953518-01 
2. 28266508— 01 
2.2 8579498-01 
2.28892488-01 
2.29205478-01 
2. 29 5 18468-0 l 
2.29811448-01 
2.3014443^-01 
2.30457428-01 
2. 307704ir~oi 
2 • ?l 083408-0 1 
2. 31 3 963 56—01 
2.31 7093 8 r -0l 
2. 120223 78-01 
2.323353o r -01 
2. 22648758-01 
2.3296133^-01 
2. 132 7432F-01 
2.33587318-01 
2 • 339C0 30 c -0l 
2. 34213298-01 
2. 34526288-0 i 
3.34635276-01 
2.35152268-01 
2.354o525*-0l 
2 • 357782 4£ -01 
2 . 3609 1 2 3 c -0i 
2.36404218-0! 
2.36717208-01 
2.37030198-01 
2.37 3431 66-01 
2.3 7* 5M 7F-01 
2.3796916F-01 
2 . 38282 15 c -0i 
2.3859514F-01 
2.3090813^-01 
2.39221 1 2F-01 
2. 395 34106-0 1 
2.3984 7098-01 
2.40160088-01 
2 • 404 7307 = -0l 
2 .40786068-01 
2 • 41 O99O5F-01 
2.41412048-01 
2.41 725038-01 
2.42038025-01 


1.491 79 10^-04 
1. 26312818-04 
i.04547L0 c -04 
8.49436568-05 
6.86118148-05 
5.62078958-05 
4.82322898-05 
4.4844 72 1 e — 05 
4.59368738-05 
5. 1125335P-05 
5.98769728-05 
7.1 350872^—05 
8.469841 5 E-05 
9. 90256728-05 
1.132880*6-04 
1.26679338-04 
1 .38324758-04 
l. 47612058-04 
1.54072535-04 
1.5 7373168-04 
1.57406645-04 
1.5427273^-04 
1.4819443^-04 
1.3951667^-04 
1.28753968-04 
1. 16 381 3Q c -04 
1.03021216-04 
8.9191266^-05 
7.53731645-05 
6.2032922 E-05 
4.94 794438-05 
3.T956079F-05 
2.7591600^-05 
1.8401536F-05 
1. 02950768-05 
3.1022257^-06 
-3. 3724850^06 
-9. 433467 i c - 06 
- 1 • 53395328— C5 
-2. 1349101 -—05 
-2. 7764954F-05 
-3 .47525048-05 
-4.2460 76 7 c -05 
-5. 1003302F-05 
-6.0392 50? C -C5 
-7.05 784885—05 
-8. I451666 C — 05 
-9 • 28 753? 08 -05 
-1.04660928-04 
- 1 • 165 7039 e -04 
-1.28389348-04 
-1.35878118*04 
-1.50856058-04 
-1.6110024^-04 
-1. 70521998-04 
- 1 ♦ 79001 5 3 8-04 
-l.8649335 c -04 
-1.9298014E-04 
- 1. 9851550F-C4 
-2.03172108-04 
- 2. C7 079 508- 04 
-2.10359878-04 
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2.42351016-01 
2.4266399=-01 
2.42976986-01 
2.43289976-01 
2 .$3602966-01 
2.43915956-01 
2.44228946-01 
2.44541936-01 
?.4485492'-01 
2.45167916-01 
2.45480906-01 
2.45 79?896-0l 
2.46] 0687F-01 
2.4641 586 c -01 
2.46732856-01 
2.4704584 c -0l 
2.47358836-01 
2.47671826-01 
2.47984816-01 
2.48297806-01 
2 .4861 079 c -01 
2 .4892378 c — 01 
2 .492 3676 c — 01 
2.4954975E-01 
2.49862746-01 
2.50175736-01 
2.50488726-01 
2.50801 71 c -01 
2.511! 4706-01 
2.5142769--01 
2.51 74C686-0I 
2.52053676-01 
2.52366666-01 
2.5267964F-01 
2.5299263 c -01 
2.53305626-01 
2.53618616-01 
2.53931606-01 
2.5424459^-01 
2.54557586-01 
2. 54e 705 76-01 
2.55183566-01 
2.55496556-01 
2.55809536-01 
’.5612252E-01 
2.56435516-01 
2.56 7485 06-01 
2.57061496-01 
2.5 73 744 8 c -0 1 
2.57687476-01 
2.58000466-01 
2.58313456-01 
2.58626446-01 
2.58939436-01 
2.59252416-01 
2.59565406-01 
2.59878396-01 
2.60191386-01 
2.60504376-01 
2.60817366-01 
2.61130356-01 
2.61443346-01 


-2.131527 36-04 
-2 .1555874 c -04 
-2.1782228E-04 
-2.19933536-04 
-2.22010446— C4 
-2.24063826-04 
-2.26095266-04 
-2.2805480F— 04 
-2.29840056-04 
-2.31329656-04 
-2.32319966— C4 
-2.32609056-04 
-2.32016546-04 
-2.30341196-04 
-2.27408966-04 
-2.23039356-04 
-2.17139846-04 
-2.09*71276-04 
-2.CC683576-04 
-1.90303706-04 
-1.78 744696-04 
-1.66313916-04 
-1.53376986-04 
-1.40470486-04 
-1.28114516-04 
-1.16917456-04 

- 1.07448 12P-04 
-1.00416176-04 
-o. 63731636-05 
-9.58207 77 c — 05 
-9 .9251758 F-06 

- 1.C6950 296-04 

-1.19O7665E-04 
-1.35615666-04 
-1.56358156-04 
-1 .80870466-04 
-2.C0597 16F-C4 
-2.38635706-04 
-2.70052646-04 
-3.01769356-04 
-3.32472106-04 
-3.60871 15F-04 
-3.65746956-04 
-4.05719726-04 
-4.19 707976-04 
-4.26769386-04 
-4.26055586-04 
-4.17105976-04 
-3.59709786-04 
-3.73993326-04 
-3.40383116-04 
-2.99604966-04 
-2.52820746-04 
-2. 01 286-04 

-1.464 .206-04 
-9.02030416-0* 

-3.40419516-05 

2.00805096-05 

7.06864706-0*’ 

1.16138246-04 

1.55247786-04 

1.87077186-04 
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2.61 7563 36— 0 1 
2. 6206932 F -01 
2.6238230*-01 
2.1:2 69529E-01 
2*630082 8F-01 
2.633212 76-01 
2*63634266-01 
2.6394I25 c -0l 
2*64260246-01 
2*645/3236-01 
2*64886226-01 
2*65199216-01 
2.65512206-01 
2.65825186-01 
2*66138176-01 
2*66451166-01 
2*66 7641 56— 0 1 
2. 6 7C 77146-01 
2*t7?90136-01 
?.67703l2 c -0l 
2.68016116-01 
2*6832510—01 
2*68642096-01 
2 *689550 76-01 
?. *9268066-01 
2.^581056-01 
2*698940-* r -01 
2. 702070^ F— 01 
2. 70*2002^-01 
2 • 708 3 3916—0 1 
2« 71 I4600 c -0t 
2* 71456596-01 
2.71771986-91 
2 . 72 0 84 976-01 
2. 72397956-01 
2. 72 7 1 G94 c -0 1 
2. 730?3^F-oi 
2*73? 36926-01 
2* 7364991 c —0 1 
2.73962906-01 
2.742 7589 c -0 1 
2. 7*588886-01 
2. 74901876-01 
2 . 752 14866-01 
2. 75527* *£-01 
2. 75t>*083 r -01 
2 • 76 1 5382 r — 0 1 
2 • 764668 l - - 01 
2. 76775806-01 
2.77092796-01 
2 .77408 78f-0l 
2.77718776-01 
2. 78031 76 c -0 1 
2. 78? 44 756—01 
2 • 7865 774--01 
2 * 789 7C726— 01 
2 « 7928371 £-0 1 
2. 78596 706-01 
2.79909696-01 
?. 8022268 — 01 
2.80535676-01 
2*80848666-01 


2 • 1 C91656 c -C4 
2.26499816-04 
2.33934236-04 
2.336491 86-04 
2*26341406-04 
2.13156196-04 
1.5522 5 74 c -04 
1 • 7409685 c — 04 
1*51 168C8 6—04 
l .27833516- C4 
1.05608536-04 
8. 560335 7 C - 05 
6.89325996-05 
5.6291788^05 
4. 81494156-05 
4.46963586-05 
4.584158C6-05 
5.116303 76-C5 
6.01232766-05 
7. 18436816-05 
8.546192 8 6-05 
1* 000616 76-04 
1.14580116-04 
1.28 I9276 c — 04 
1.400091 3E-04 
1.4541601 c — 04 
l .5553 391 c — 04 
I.5924il4£-C4 
l • 592442 1 6-04 
1.56041676-04 
1 • 498716 l 6—04 
1.41100156-04 
1.30260176-04 
1.1 783663 c -04 
1.04468606-04 
9.06 77C 7 7 p -C c 
7.6931 l 656— C5 
6.3714066F-C5 
5. 13186286-05 
3 .89866 88 6—05 
2.5826 144 c -05 
2.C8453 I 5E-C5 
1.29189066-C5 
5 .881962 8E-C6 
-4. 64497^3 6— C 7 
-&.4575551E-CF 
-1.2*235966-0* 
- 1.85265136-C* 
-2.51057406-05 

- 3. 2352 101 £-05 
-4.C3666026-C5 

- 4 . 92946 7 6 £-05 
-5.51 1426 1 r -0 c 
-6.S77982H6-05 

- 8. 115575* c -o c 
-9. 30 725096-05 
-1.05351 7 2 £— C 4 
-1.1 f 735486-04 
- 1 . 2 ^ 97852^-04 
-1.41865826-04 
-1.53148266-04 

- 1.63640226-04 
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7.81161656-01 
2.8 147464E-01 
2.81 787636—01 
2. 8210061 E-01 
2. 8241360? -01 
2.8 2 726 5 96 -01 
2.8303958 c -01 
2.83352S7E-01 
2 • 83664568—01 
2.83978556-Oi 
2.04291546-01 
2.8460453^-01 
2.8491 752?— 01 
2.8523051 6—01 
2. 855434^6-91 
2.85856486-01 
2.86169476-01 
2.06482466-01 
7. 86 7554 5 P-0 1 
2.87? 30448-01 
2.87421438-01 
2 .07734426-01 
2 . 88047415-01 
2.80360406-01 
2. 086 73386-01 
2.885863 /P-01 
?. 89295368-01 
2 • 89M2 15 C -01 
2. 09925346-01 
2 .902 303 3 6—01 
7. 9055 1 326-01 
2 . 9C0643 1E-0 1 
2.511 / 73 0F-01 
2 .91490? 9E-01 
2.9I0O320 C -O1 
2.921 1626E-01 
2.52429756-01 
2.52742246-01 
2.93055736-01 
2 .933682? c -01 
2.93681218-01 
2.93594206-01 
2.9^30/196-01 
2.546 20186-01 
2.945331 7*-01 
2.9S24615E-01 
2 • 55559146—0 l 
2.95872136-01 
2.56105126-01 
2.9649811E-01 
2.96811106-01 
2.97124096-01 
2 .974 3 7086—01 
2. 977500 7*-0i 
2.98063066-01 
2.98376056-01 
2.586 890 36-01 
2 .99002026-01 
2 * 993 1501 r -OI 
2. 99628006-01 
2.9994C99 C -01 
3.00253986-01 


-1.73229456-04 

- 1.81797746—04 
-1.89278356-04 
-1. 95682755-04 
-2.01061 C86— 04 
-2.05488346-04 

- 2.09091 8Q C — 04 
-2.11994796-04 
-2.14378276-04 
- 2. 16 3 8182 E— 04 
-2.18151166-04 
-2.19804355-04 
-2.21431 586-04 
-2.2 3081 10 C -C4 
- ’.24760926-04 
-2.264296 35-04 
-2.2800 1876-04 
-2.29350586-04 
-2.3030933—04 
-2.30718926-04 

- 2. 303286 35— 04 
-2.28553976-04 
-2.26405126-04 
-2.7253 79 /E-C4 
-2.17181216-04 
-2.1032 7446-04 
-2.0l97950 c -04 

- l.9222602 c -C4 

- 1 .8126 74C 6-04 
-1.69383606-04 
-1. 56935556- C4 

- 1. 4436521 6— 04 
-1.32250055-04 
-1.21 1 5240 E- 04 
-1.1 I67826 c -04 
-1.04416386-04 
- 1.00045 12 P—04 
-9. 50741816-05 
-1.01956826-04 
-1.0902042^-04 
-1. 20485 006-04 
-1.36354335-04 
-l. 56478046-04 

- 1. 80447 34 c - 04 

- 2. 076 9 72 5 p-04 
-2.3750114E— 04 
-2.680O9/0 C -O4 
-3.00543596-04 
-3.31 5 5449 P— 04 
-3.60395546-04 

- 3.0 5818 ICE— 04 
-4.0661530E-04 
-4.2146831 F— 04 
-4.29424086-04 
-4.79690716-04 
-4.2164648P-04 
— 4. C 50 75116- 04 

- 3.80047816-04 
-3.46960 I 06-04 
-3.06486526-04 
-2 .59 77734* -04 

08092506-04 
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3.00566975-01 
3.00879965-01 
3.01192958-01 
3 • 01 5 05 94 c -0 1 
3.01 818928-01 
3.021 3191 c -01 
3 .02444905-01 
3.02757895-01 
3.030 7C 888-01 
3.03383878*01 
3.03696868-01 
3 .04009 8 5 F— 01 
3.04322849-01 
3.04635838-01 
3. C4948828— 0 1 
3.05261808-01 
3.05574798-91 
3 . 05 88 778 c -01 
3.0b?0077F-01 
3.0651376^-01 
3.06826755-01 
3.07139748-01 
3.0/452738-01 
3.0776572^—01 
3.08078 71 r -0i 
3.03391698-01 
3.0«70468 c -01 
3.090!767r-0! 
3.0933366^-01 
? •0964‘ a 65F— 01 
3.09956648-01 
102 6 96 3 8-01 
3.105*2628-01 
3 • 1089561^-01 
3.11208608-01 
3.11521585-01 
?. 1.034578-01 
3.121 4 756 P — 01 
3 . 12460558-01 
? .1277354F-01 
3.13086535-01 
3.13399528-01 
3. 13712515-01 
3.14025505-01 
3.1433049P-OI 
3.1465148^-01 
3 . 14964468—01 
3. 1527745 — 01 
3.15590448-01 
3.15903438-01 
3.16216428-01 
^.16529418-01 
3. 16 84240C- 01 
3 • 1 71 5539F-0 1 
3.1 74683 8 C — 0 1 
3.17781375-01 
3.1809435»— 01 
3.18417348-01 
3.18720335-01 
3 • 1903332 c -01 
3.19346315-01 
3 . 1965930^-01 


- 1 • 52 869 5 5 F- 04 
-9.59782695-05 
-3.90018615-05 
1.61074025-05 
6.77242215-05 
1 .14 384 8 2 p -04 
1.54636345-04 
l • 8754847 c -04 
2.12382295-04 
2.28843738-04 
2.36952078-04 
2.3 7164498-04 
2 .301492 7 c -04 
2.17055225-04 
1.99102725-04 
1.77646625-04 
1.54335098-04 
! .30517955-04 
1.0 7755715-04 
8.72061355-05 
7.001 22535-C5 
5. 690 74 7 fc* 1 - 05 
4.83813165-05 
4. 46453065-05 
4.56168635-05 
5.08747608-05 
5.9863851 c -05 
7.17 CQ9055— 05 
~ * 55307378—05 
1.00380225-04 
1. 1516938F-04 
1. 29059505-04 
1.41 140468-04 
1.507826 7*- 04 
1.57498258-04 
1 .609600 38-04 
1.61074408-04 
1.5 7949255-04 
1.51800745-04 
1.43062335-04 
1. 32231328-04 
1.19819425-04 
1. 06482528-04 
9.2 7044465 - C5 
7.9C5370 7 F-05 
6.5935522 E-05 
5.36654658-05 
4.?486053 c -05 
3.24858865-05 
2.3659735^-05 
1.58851575-05 
8.99642935-06 
2.70804075-06 
-3.25907728-06 
-9.1 7378938-06 
-1.53710745-05 
-2.20796755-05 
-2.95016495-05 
-3.781 I4405-C5 
-4. 70644635-05 
-5.72621115-05 
-6.83435145-05 
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3.1997229E-01 
3. 20285286-01 
3.20558276-01 
3.2091126E-01 
3.2122425F-01 
3.21537236-01 
3.2185022^-01 
3.2216321F-01 
3.22476206— 01 
3.22789196-01 
3.2310218E-01 
3.2341517F-01 
3.2372816E-01 
3.24041 l 5E-01 
3.2435414E-01 
3. 2466 71 2 E— 01 
3.249831 IE-01 
3. 252931 0E-01 
3.2560609'-0l 
3. 2591 908E-01 
3.2623207E— 01 
3.26545066-01 
3 .2685805 c — 91 
3.271 7104E-01 
3.2748403F-01 
3. 2779702^-01 
3. 281 1000F-01 
3.2342299F-01 
3 .2873598^-01 
3.2904897E-01 
3. 2936 1 96 E— 01 
3 .296 7495F-01 
3.2998794P-01 
3. 3030093E-01 
3. 3361 392 *—01 
3. 3092691E-01 
3.312 3989F— 01 
3.3155288F-01 
3. 318658 7E-01 
3.321 7886 r -01 
3 • 32491 85E-01 
3. 3280484F-01 
3 .331 1 783E-01 
3.3343082E-01 
3 . 33 74381E-01 
3 • 3405680E-01 
3.^436979F-01 
3.3468277E-01 
3. 34995T6 c -01 
3.35308756-01 

3.3562174F-01 

3.3593473«=-01 
3.3624772F-01 
3.3656071 c -0l 
3 . 36873 70E-01 
3 • 371 8669F-0 l 
3.3 749568E-01 
3.3781266F-01 
3.3812565E-01 
3 • 3843 864F-01 
3. 3875163 c -01 
3 • 3906462E-01 


-8. 0174203E-05 
-9. 25793 1 7F-05 
- 1.0532584E-04 
- 1.181 71 05F-04 
- 1 . 30855146-04 
-1.43144536-04 
-l. 547795 8F-04 
-1.65584686-04 
- 1 . 7539902 c -04 
-1.84104216-04 
-1.9165266F-04 
-1.58041606-04 
-2.03311266-04 
-2.07555856-04 
-2. 1090596 E-04 
-2. 13501 8 7E- 04 
-2. 1550811 E-04 
-2.17093286-04 
-2.1840655F-04 
-2.19602386-04 
-2.2C78104F-C4 
-2.22006876-04 
-2. 2330744 F-04 
-2. 24662 16 E-04 
-2.2601034F-04 
-2.2 7220 7CF-04 
-2.2816571 E-04 
-2.2664585F-04 
-2. 28469436-04 
-2.2744705F-04 
-2.2 53 25 15 E-04 
-2.2198571 E-04 
-2. 172 7109 c - 04 
-2.11 11131 E-04 
-?.C347826 c -04 

- 1.5445034E-04 
- 1.8418308F-04 
-1.7289753F-04 
-1.6098564f-04 
- 1 .4885 1 99F-04 
-1.3699619E-04 
-1.2594625F-04 
-1. 16381836-04 
-1.08895586-04 
-1.040433 7E-04 
-1.02 46951 6-04 
-1.0461635^-04 
-1 .108517CE-04 

- 1.21412786-04 

- 1 . 36 38 197 f-04 
-l. 55624496-04 
-1. 7885531F-04 
-2.05467266-04 
-2.34 76 77 1 F-04 
-2 .65893656-04 

- 2.9763 759F- 04 
-3.2884203F-04 
-3.58278556-04 
-3.844 754*6-04 
-4.06199096-04 
-4.22285226-04 
-4.31526786-04 
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A* UJ IhJ 


3* 3937761 £—01 
3.39690606-01 
3.4000359E-91 
3*403 1658 e -01 
3.40629576—01 
3.4094256*— 01 
3.41255546— 01 
3*41566536-01 
3.41661526-01 
3.4219451^-01 
3. 425075 06-01 
3.42 820496-01 
3.43133486-01 
3.4344647^-01 
3*43755466-01 
3.440 72*, 56-01 
3 • 443 8543 p — 01 
3 . 4469842^-01 
3 * 450 1141 6—0 1 
3. 4532440 r — 01 
3.45637396-01 
3.45950^8^-0 l 
3.46263376-01 
3. 46 5^6 366—01 
3.46685355-01 
3 • 472 0234F-01 
3.4751533^-01 
3.47828316—01 
3.481413 0P-01 
3. 48454296-01 
3*48767286-01 
3.490802 7 C -01 
3.49393266-01 
3 .49706256-01 
.50019246-01 
.503^2236-01 
•5064522E-01 
.50958206-01 
3.512 711 9 C — 0 1 
3.51584186-01 
3.51 8971 7F-0 1 
3.5221016E-01 
3. 52523156-01 
3.52836146-01 
53149136-01 
3 .5346212^-01 
3.53775116-01 
3.5408810^-01 
3.54401186-01 
3.54714076-01 
3.55027066-01 
3.55340056-01 
3. 5555304 e-oi 
3.55966036-01 
3.56279026-01 
3.56592016-01 
3 .5690500^-01 
3.57217996-01 
3.57530976-01 
3.57843966-01 
3.58156956-01 
3 . 58469946—01 


-4.33121796-04 
-4.26470606-04 
-4.11201176-04 
-3.07318936-04 
-3.55200586-04 
-3*1 550 736 c -04 
-2.69167306-04 
-2. 17595636-04 
-1.62170746-04 
-1.04747486-04 
-4.70485456-05 
9.21439066-06 
6.21049576-05 
1.1018215E-04 
1.51934906-04 
1.86328516-04 
2. 1262765^-04 
2.30402326-04 
2. 39643576—04 
2.40 768?b c -04 
2 . 3447021 c - 04 
2.21726556-04 
2.03927736-04 
1.82370606-04 
1.56751946-04 
1.34450046-04 
1.11087836-04 
0.99OOfc5O c -O5 
7. 194965? F— C5 
5.01543756-05 
4. 89920655-05 
4.47278016-05 
4.52839926-05 
5. 02604*1 *-05 
5.99661016-05 

7 . r '«* 10 586-05 

8. <.*,©61646-05 
9. 99280926— C5 
1.15001926-04 
1.29224626-04 
1 .4166 725 e -04 
1.5L670306-C4 
1 • 5073061 c -04 
1.6250396E-04 
1.62891336-04 
1 *5998398 c -04 
1.54033206-04 
1 .4541297 c -04 
1.34670926-04 
1.22339016-04 
1.09053706-04 
9.53298466-05 
8.1736 745 c -05 
6.86882236-05 
5.65095176-05 
4.54278406-05 
3.55305276-05 
2.60211666-05 
1.41631206-05 
1.23864396-05 
6 . 18281116-06 
2.74282496 -u7 
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3 • 58782938—01 
3.59095926-01 
3. 594O091F— 01 
3. 59721908-01 
3 .60034 89 c -01 
3.60347888-01 
3 .60660 87 £-01 
3*60973856—01 
3. 61286046— 01 
3.61555838-01 
3. 61 91 2828-01 
3 . 6222581^-01 
3.62538808-01 
3.62051796-01 
3 .631 6478E-01 
3 • 634 7777F— 01 
3*63 7907oE-0l 
3.641 0374E-01 
3.64416736-01 
3 • 64 72972E -01 
3.6504271 c — 01 
3.65355708-01 
3. 65668698-01 
3. 65981688-01 
3.6629467E-01 
3 . 6660 ?66 c -0 1 
3.6692C65E-01 
3. 672 33646-01 
3.67546628-01 
3.67855618-01 
3.681 7260P-01 
3.68485598-01 
3 .6879858^-01 
3.69111578-01 
3. 694245 6P -01 
3. 69737558-01 
3.70050548-01 
3. 7036353E-01 
3 . 70676516-01 
3. 70989508-01 
3.71302498-01 
3. 71615486-01 
3.71928478-0! 
3.7224146E-01 
3.7255445P-01 
3* 72 8t»7448-01 
3 . 731 8043 c -01 
3.73493428-01 
3 • 7380640E— 01 
3.74115398-01 
3.74432388-01 
3. 74 7453 7^-01 
3 • 75058368-01 
3.75371358-01 
3.75684348-0! 
3. 75997338-01 
2.76310328-01 
3.76623316-01 
3.76936306-01 
3.77249286-01 
3.77562276-01 
3.77875266-01 


-5.62301376-06 
-1.1851940E-05 
-1.86 4598 58—05 
- 2. 62 1 33 5 7 F— 05 
-3.472 7558F-G5 
-4.42443888-05 
-5.4/702556-05 
-6.62195776-05 
-7.84631156-05 
-9.13001228-05 
- 1 .0450021 p — 04 
-1.1 780304 p- 04 
-1.30955378-04 
-1.43669636-04 
-1.55701206-04 
-1.66848736-04 
-1.76945486-04 
-1.85870536-04 
-1. 935643 8E-04 
-2. 00028006-04 
-2. 05279 3 5 s - 04 
-2.09419656-04 
-2.125785CP-04 
-2.14894998-04 
-2.16569706-04 
-2.17767056-06 
-2.13660036-04 
-2.1940267 c -04 
-2. 20 12364E-04 
-2.20911938-04 
-2.21807238-06 
-2. 22820426-04 
-2.23894686-04 
-2.24954668-04 
-2.25864326-04 
-2.26422278-04 
-2.26443888-04 
-2.25 740848-0- 
-2.24097076-04 
-2. 21354286-04 
-2.17324396-04 
-2.11908666-04 
-2.05069346-04 
-1.96845386-04 
-1.97347206-04 
-1.76 798876-04 
-1. 65494 71P-04 
-1.53790138-04 
- 1. 422 37 79F-04 
-1.31331406-04 
-1.21621156-04 
-1.1 3 79855 F-04 
-1.08448186-04 
-1.06135356-04 
-1.07346656-04 
-1.12551416-04 
-1. 2 20 1 36 P 6-04 
-1.358356, 6-04 
-1.53967356-04 
-1.76141146-04 
-2.01924146-04 
-2.30579006-04 
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3 • 781 0825F-O l 
3 •78501245-01 
3.78814236-01 
3 • 7912 722--01 
3. 79440216-01 
3 • 7975320 c -0l 
3.80066196-01 
3.80379176-01 
3. 80692166—01 
3.81005156-01 
3.01310146-01 
3.61631136-01 
3.8194412^-01 
3 . 822571 1F-01 
3. 82570106-01 
3 .82883096—01 
3.03196006-01 
3 • 83509016-0 1 
3. 83022056-01 
3.84135046-01 
3.84448036-01 
3.04761026-01 
3.85074016-01 
3.85387006-01 
3.85699996-01 
3.86012986-01 
3.86325976-01 
..86638966-01 
3. 86951946-01 
3.87264936-01 
3.87577926-01 
3.8789091^-01 
3.80203906-01 
3 • 8851 6896-0 i 
3 .88929886— 01 
3.89142876-01 
3.8945586F-01 
3.89768056-01 
3.90081846-01 
3.90394826-01 
3 .9070701 c -0 1 
3.91020006-01 
3.91 33 3 796-0 L 
3.91646786-01 
3.91959776-01 
3.922 72 766-01 
3.92585 756-01 
3.9209874^-01 
3.93211736-01 
3.93524716-01 
3.53837706-01 
3.94150696-01 
3.94463686-01 
3.94776676-01 
3.95089666-0 1 
3 .954026 5F-01 
3. 95715646-01 
3.96028636-01 
3* 96? 41 626-01 
3.966 54616-01 
3.96567596-01 
3.5728058^-01 


-2.61253736-04 
- 2.92972526-04 
-3. 24419H 6—04 
-3.54355716-04 
-3.81514186-04 
-4.04421216-04 
-4.21886286-C4 
-4.32825716-04 
-4.36187156-04 
-4.31 33068°— 04 
-4.178589 76-04 
-3.55660266-04 
-3.64963 756-04 
-3.26447976-04 
810183 4F- 04 
- 2.29831 736—04 
-1.74542726-04 
-1.16662316-04 
-5.81 056006—05 
-0.02879656-07 
5.35412016-05 
1.03243916-04 
i. 46894736-04 
1.832C3096-04 
2.11 38553E-34 
2.30965986-04 
2. 41852056-04 
2.44320176-04 
2. 39111^76-04 
2.2 712945E-04 
2.09762256-04 
1.88367736-04 
1.64490146-04 
1.39811406-04 
1.15714176-04 
9.37104576-05 
7.48543866-05 
6.01249546-05 
5.00785506-05 
4. 500 8552 05 

4.4863446^-05 
4.53372016-05 
5.77935736-05 
6.55066666-05 
8.3434558F-05 
9*8651351 E-05 
l .14014036-04 
1.28586606-04 
l * 4 1 52301 F-04 
1. 520 1574E-04 
1 . >9 5 7 ? 1 9 E - 0 4 
1.63^26306-04 
l .64646 87F-04 
1.62119246-04 
1. 56408046-04 
1.48119576-04 
1.3757657^-04 
l .25376296— 04 
1.12104106-04 
9. 85483 126-05 
8*45 7652 5F-05 
7.19736 72 c -05 
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3.97593S76-01 
3.97906565-01 
?.982 19556-01 
3.58532546-01 
3. 9884553^-01 
3. 99156526-01 
3.994715 l c — 01 
3.99764506-01 
4*00097486—01 
4 *004104 76—01 
4.00723466-01 
4.0103645^-01 
4.01349446-01 
4.0166243*— 0 I 
4.C197542F-01 
4.02288416-01 
4.026 01406-01 
4.02914396-01 
4.03227386-01 
4.03540366-01 
4.0385335^-01 
4.04166346-01 
4.04479336-01 
4.04792326-01 
4.05105316-01 
4. 05M 8306-01 
4.05731296-01 
4.06044286-01 
4.06357276-01 
4.066 7°256— 01 
4. 065 83 2 46— 0 1 
4.07296236-01 
4.0760922—01 
4. C7922216-01 
4.082352 06-01 
4.C8548196-01 
4.38861 186—01 
4.091 74176-01 
4.09487166-01 
4 .098001 5 C — 01 
4.101 131 36-01 
4. 1042612 s — 01 
4.10739116-01 
4.11052106-01 
4 • 1 1 36 509 6—01 
4.11678086-01 
4. 1109 1 076- 01 
4.12304066-01 
4.12617056-01 
4.12930046-01 
4.13243026-01 
4.13556016-01 
4.13669006-01 
4.14181996-01 
4.14494986-01 
4. 1480797 c — 0 1 
4.15120966* *1 
4.15433956 l 
4.15746946 01 
4.16059936-01 
4.16 3 72926-01 
4.16685906-01 


5.98266376-05 
4.87919686-C5 
3.89543756-05 
3.031 771 IE-C5 
2.27468466-05 
1.60487126-05 
9.92573236-04 
4.07901046-06 
-1.77437786-06 
-7.99023896-06 
-1.48159366-05 
-2.24622811-05 
-3.11088176-05 
-4.00149506-05 
-5.15820756-05 
-6. 33323736-05 
•592781 46— C5 

- 8.9 16231 56— C 5 
-1.02 794556-04 
-1.16559106-04 
-1.30178966-04 
- 1.43 3443 66- C4 
—I. 5581 582E-04 
-1.67367296-04 
-1.77839406-04 
-1.87076216-04 
-1.95004CI6-04 
-2.0161690^-04 
-2.069242*6-04 
-2. 11045226— C4 

- 2. 140967 26—04 
-2.16229136-04 
-2.1 7625 70 c — 04 
-2.18475556-04 
-2.1895708^-04 
-2.19261586-04 

- 2. 195189J5—04 

- 2. 198485 56-04 
-2.20318906-04 
-2.20954946-04 
-2.2174001 E-04 
-2.2259643^-04 
-2. 233903 96-C4 
-2.23989906-C4 
-2.24205086— 04 

- 2.23848346— 04 
-2.22687226-04 
-2.20529386-04 
-2.1 7206MF-04 
-2.12611916-04 
-2.06641476-04 
-1.99297526-04 
-1.90677106-04 
- 1 . 809 1921 c -04 
-1.70309346-04 
-1.59190706-04 
-1.4 7993 156-04 

- 1.37190146-04 
-1.27421976-04 

- 1 • 1 92 5 7006—04 
-1.132673 76-04 
-1.10127476-04 
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4.169S889=-01 
4.17311886-01 
4.17624876-01 
4.1 793786E— 01 
4.18250856—01 
4 -18563846-01 
4.i8»76B3F-0J 
4.19189B26 ' I 
4.1450281 6—01 
4.19815796-01 
4.20128786-01 
4.20441776-01 
4. 20 734 ?6 C -01 
4 .21067751-01 
4.21380746-01 
4.2l643~36-0i 
4.22006726-01 
4.22319716-01 
4.22632706-01 
4.22945696-01 
4.73258676-01 
4.75571666-01 
4*23884656-01 
4.24197646-01 
4.245106 3 6-01 
4.24823626-01 
4.2513661=— 01 
4.75440606-01 
4.75762596-01 
4.26075586-01 
4.26388566-01 
4.24701556-01 
4.27014546-01 
4.273275 3 r -01 
4.27640526-01 
4.27953516-01 
4.28266506-01 
4.285 79496-01 
4.28892486-01 
4.29205476-01 
4.2951 o46 r -01 
4.29831446-01 
4.3014443?-01 
4.30457426-01 
4.30770416-01 
4.31083406-01 
4.3139639F— 01 
4.31 70938 c -01 
4.32022376-01 
4.323353*6-01 
4.32648356-01 
4.32961336-01 
4 . 332 7432"!-01 
4.335873 1 F -01 
4.3390Q?0 C -01 
4. 34713^96-01 
4.34526286-01 
4.34839276-01 
4. 35152266-01 
4.35465256-01 
4.35778246-01 
4.36091236-01 


-1.10323526-04 

-1.14314096-04 

-1.22419766-04 

-1.34859226-04 

-1.51615316-04 

-l.7252l50 c -04 

-1.97163846-04 

-2.2492985=-04 

-2.55070C06-04 

-2.86481756-04 

— 3. 180 766 IE— 04 
-3.48632126-04 
-3. 7667346 =-04 
-4. 00989576-04 
-4. 20141966-04 
-4.3298719=— 04 
-4.38561526-04 
-4.35982756-04 
-4.24786206-04 
-4.04 764 726-04 
-3.76093996-04 
-3.39222926-04 
-2.55094546-04 
-2.44837196-04 
-1.89821966-04 
-1.31 99079F— 04 

- 7.266888 1=-05 
-1.42250866-05 

4.16757426-05 
9.33618266-05 
1. 3912321 c -04 
1.77803626-04 
2.C834982E— 04 
2.30231 70 c -04 
2.432 50=0=— 04 
2.4764561=- 04 
2. 43955746-04 
2. 3322287=- 04 
2. 1 6 554? 66-04 
1.955452 36-C4 
1.7161870 c -04 
1.46546436-04 
1 »2I8206fi=-04 
9.88229596-05 
7.89238786-05 
6 .30000416-05 
5. 1 795848=— 05 
4. 5624 762 =-05 
4.45064716-05 
4. 819543 7=-05 
5.6C73016E-C5 
6.7 3820 3 7=— 05 
8.12267956-05 
9.64737146-05 
1.12153I2 c -04 
1.2715618E-04 
1 .4063368=— 04 
1.51744626-04 
1.59954036-04 
1 .64851046— 04 
1.66293656-04 
1.64 306 51 =-04 



4. 3640421 6 -Ot 
4.367i720 c -0l 
4.37030196-01 
4. 373431 86-01 
4. 37656176-01 
4.3 7969166—01 
4.3828215 c -01 
4.38595146-01 
4.38908136-01 
4. 39221 1 2E— 01 
4. 39534106— 01 
4. 3984 709 c — 01 
4.40160086-01 
4 .404 73G7 C — 01 
4.40786066-01 
4.41099056-01 
4.41412046-01 
4.417250*6-01 
4.42038026-01 
4.42351 016—01 
4.4266399 c -01 
4.42976986-01 
4.43289976-01 
4. 4360296 c — 01 
4.43915956-01 
4.44228946-01 
4.44541936-01 
4.4485492 c -0l 
4.45167916-01 
4. 4548090^-0 | 
4 .45 793896—01 
4.46106876-01 
4.464 1 986-t-O 1 
4.46 732856-01 
4.47045846-01 
4.47359836-01 
4.47671826-01 
4.4/98401E-OI 
4.48? 97806—01 
4 .486 10 796-01 
4.48923786-01 
4.49236766-01 
4.49549756-01 
4.49862746-01 
4.50175736-01 
4.50488726-01 
4.50801 71 c -0 1 
4.51114706-01 
4.51427696-01 
4.5174C686-01 
4. 52053676-01 
4.52366666-01 
4.52679646-01 
4.52992636-01 
4.53305626-01 
4. 5361 86i c -01 
4.53531606-01 
4 . 542 44 59 F- 01 
4.54557586-01 
4.54870576-01 
4 . 551 8356 c -01 
4.55496556-01 


1.59154456-04 
1.51 171846—04 
l. 40919606-04 
1.28962196-04 
1.15874306-04 
1.02 324396-04 
8.87837646-05 
7.57763036-05 
6.36296736-05 
5.25788666-05 
4.27523866-05 
3.41365906-05 
2.660016? 6—05 
1.994324 76-05 
1.39122016-05 
8. 1543194—06 
2.37532906-06 
-3.77641456-06 
-1.C5 728476-05 
-1.82224546-05 
- 2.692004 7C-05 
-3.67332226-05 
-4.7663972^-05 
-5.96449556-05 
-7. 2527310^-05 
-8.61 233 39 F- 05 
-1.00155106-04 
-1. 1435301^-04 
-1.20439436-04 
- 1 .42092 75F-04 
- 1 .55053696-04 
- 1.670861C6-04 
-1. 78C003 C 6-04 
-1.0 76182 36-04 
-1.55877706-04 
-2.02749436-04 
- 2. 00246 4 O c — 04 
-2.12450276-04 
-2.15481916-04 
-2.175000 7P-04 
-2.1867670^-04 
-2.1924508^-04 
-2.19362396-04 
-2.19245126-04 
-2.19047816-04 
-2. 10907626-04 
-2.18924116-04 
-2.19135316-04 
-2.1956834C-Q4 
-2.20156896-04 
-2.2081451 F-04 
- 2 * 2 1 4242 46- C4 
-2. 21777956-04 
-2.21690546-04 
-2.20975296-04 
-2. 19440096-04 
-2.16858316-04 
-2.1 33951 5P-04 
-2.08047656-04 
-2.01666316-04 
-1.9399 7436-04 
-1.651 5621 6—04 
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4.55809536-01 
4.5612252E-01 
4.56435$1F*0! 
4.56 748506-01 
4.57061496-01 
4.573 74486-01 
4.57687476—01 
4. 580004 6f -01 
4.58313456—01 
4. 586 2644 E— 01 
4.58939436-01 
4.59252416-01 
4.5954540 c -0l 
4.598 78396-01 
4.6019138^-01 
4.60504376-0! 
4. 60817366-01 
4.61130356-01 
4.61443346-01 
4.61756336-01 
4 .6206932C-O1 
4. 62 3823 0 r -01 
4.6269529F-01 
4.630082 8 £-3! 
4.63321 ? 7**-01 
4.636 34766 -Cl 
4. 639472 5 C — 01 
4. 642 602 4f -01 
4. 6457? 2 3 C - 01 
4.648 8*22 r -01 
4 . 651 ^ 971^.01 
4. 6551220 r -01 
4,658251 8F -01 
4.661381 76—01 
4.66451 16E—01 
4.66764156-01 
4.67077145-01 
4. 67390136-01 
4.677031 2 r -01 
4.68016116-01 
4.68329106-01 
4.68642096-0! 
4.68955076-01 
4.69268066—01 
4.69581056-01 
4.6989404 r -0l 
4. 7020703F-01 
4.7052002E-01 
4. 708330 i F-01 
4.71146006-01 
4. 7145899 C -01 
4.71771987-01 
4.7208497^-01 

4.72 39 795F-01 
4.72710946-01 
4. 7302393 c -01 
4.73336926-01 
4 • 736499 1 c — 01 

4.73 562906-01 
4. 742 75896-01 
4. 7458888F— 01 
4. 7490! 8 7 P-0 1 


- 1.75351 16 c - 04 
-1.64857696-04 
-1.54109976-04 
-1.435631 IE— 04 
-1.33753246-04 
-1.25220796-04 
-1. 18654226-04 
-1.146! 11 3 c -04 
-1.13625556-04 
- 1. l6?5766 c -04 
-1.22869806-04 
-1.33687016-04 
-1.4881 855P- 04 
-1.681563 76-04 
-1.9141821P-04 
-2.18034776-04 
— 2.4728758P— 04 
-2.78311656-04 
- 3.09919146—04 
-3.4087352 c — 04 
-3.69949466-04 
-3.95629826-04 
-4.16639206-04 
-4.31791076-04 
-4. 39850726-04 
-4.39993166-04 
-4.31592306-04 
-4.14311586-04 
- 3. 88 16982 c -04 
-3.53597316-04 
-3.11313186-04 
-2.6233601^-04 
-2.081 8544F-04 
-1.50 320 7 7^-04 
-9.073692 8t>05 
-3.11472476-05 
2.63 7436 76-05 
8.CC640C6F-05 
1.28353006-04 
1.6969334E-C4 
2.03108336-04 
2.27811086-04 
2.4352li2 c -04 
2.50352026-04 
2.48803256-04 
2.396 7551 r -04 
2.2427695^-04 
2* C 388056 c - 04 
l • 802 30 7 76-04 
1 • 5481 30 8 p— 04 
1.2S41 2*06-04 
1.05421 3C6-04 
8 • 42 7434 2 c - 05 
6.69802336-05 
5.4 3126926— C5 
4.672961 76-05 
4.43433746-05 
4.69092635-05 
5.39662886-05 
6.46849056— C 5 
7.8229472E-05 
9.34554746-05 
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4.7521484E-01 
4.7552784F-01 
4. 7584083 ,I -0l 
4.74153826-01 
4.76466816-01 
4.76779806-01 
4.77092796-01 
4. 77405 786-01 
4.7771B77F-01 
4.78041 76 c -01 
4.78344756-01 
4. 7865 774 c — 01 
4.78970726-01 
4.7928371 c — 01 
4 . 795 96 70 e -0i 
4 . 79909fe9f — 01 
4. 8022268 c -01 
4. 805 3 56 7F -01 
4. 8084866 e -01 
4 . 81 16165F— 01 
4 . <*1474646— 01 
4.81 787636-01 
4 . 82 1 0061 F— 01 
4 . 92413606-01 
4. p? 726596-01 
4. 83039584—01 

4.8335257F-01 
4. P3665566— 01 
4 . 639 7855F-C1 
4.8429154F-01 
4.«460453 c -01 
4 . 8491 752 F -01 
4.0523051 F-01 
4.35543496-01 
4.85956486-01 
4. 661694 7 C -01 
4.86482466-01 
4.8679545 c -01 
4. 871 0844F-01 
4.87421436-01 
4. 87 7344 2 f— 01 
4.88047416-01 
4.83360406-01 
4.88673386-01 
4.885863 7C-01 
4.89299366-01 
4.8961235=-01 
4.8992534E-01 
4.90238336-01 
4.90551326-01 
4.90864316-01 
4.91 1 ?730 c -01 
4. 5L49C 296-01 
4.9180328F-01 
4.92116266-01 
4.52429256-01 
4.92 74224F-01 
4.93055236-01 
4.9336822 c -01 
4.93681216-01 
4. 93 9942 0E— 01 
4.94307196-01 


1.0931963E-C4 
1.24807886-04 
1.3886546^-04 
1.50746486-04 
1.59748566-04 
1.654826C6-04 
1.67723466— C4 
1.66475296-04 
1.61 95988F— 04 
1. 545306 16-04 
1.44669046-04 
1. 330132 76-04 
1.20120546-04 
1.06659806-04 
9.31658406-05 
8.01081216-05 
6. 7914 79 76 -05 
5.67937636-05 
4.6909610^-05 
3.6250952E-05 
3.07219446-05 
2.41012966-05 
1.81338876-05 
1.2462038F-05 
6.7684213E-06 
7.49942 83F-07 
-5.94301466-06 
-1.35379406-05 
-2.21833906-05 
-3.20079166-05 
-4.301 0680F-05 
-5.51263696-05 
-6.8222757F-05 
-8.21005916-05 
-9.64887496-05 
-I.lll050e6-04 
-1.256540*6-04 
-1.39831906-04 
-1.5331892E-04 
-1.65884506-04 
-1.77320986-04 
-1.87436926-04 
- 1.96 143 82 E- 04 
-2.0338214E— 04 
-2.091 7183 c -04 
-2.13569806-04 
-2. 16682126-04 
-2.18690916-04 
-2. 19762046-04 
-2.20105806-04 
-2.19922736-04 
-2. 1941388 e -04 
-2.13765756-04 
-2. 1815191C-04 
-2.1768630E-04 
-2. 17455126-04 
-2.17471316-04 
-2.17747616-04 
-2.1 81 8448F-04 
- 2. 18692 3 3F— 04 
-2. 19121626-04 
-2.15271646-04 
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4.94620186-01 
4.94 93317E-01 
4.552461 5F-01 
4.9555914*-01 
4.95072136-01 
4.96185126—01 
4 .96498 HE— 01 
4.96811106-01 
4.971 2409 c — 01 
-..9743708E-01 
4 .5 7 750076— 01 
4.98O6306F-01 
4.98376056-01 
4.9868903E-01 
4 . 9900?02F*0l 
4 .993 1 501 c — 01 
4.996 2 80 OF— 01 
4.5994099^-01 
5.C025398E-01 
5.00566976—01 
5 • OOtf 79966—01 
5*011 92956—01 
5.0150594^-01 
5.0101*926-01 
5.021 31 916-01 
5 • 0244490F— 01 
5.0275789P-01 
5.03O7C886-01 
5.0338387^-01 
5.03696866-01 
6.0400985F-01 
5.04322846-01 
5 .046 358 3E— 01 
5 • G494882 c — 01 
5.05261806-01 
5.0557479F-01 
5.05687786-01 
5 • 06? 007 7F— 01 
5.0651376 c -0l 
5 .06826 75E* 01 
5.071 39 74F-01 
5.0 7452 736-01 
5. 0776572F— 01 
5.C8O7071F-OI 
5 .08391 69F— 01 
5.0870468F-01 
5.09017676-01 
5 . 093 3066F— 01 
5.09643656-01 
5 *0995664 E— 01 
5. 10269636-01 
S.1058262F-01 
5. I 08956 l c — 01 
5.11 20860E-01 
5.11521 586-0 1 
5 . 11 8345 7F-01 
5.121 4756F— 01 
5 . I246055 c — 01 
5*12 77354F— 01 
5.13086536-01 
5. 13 39952E-01 
5.13712516-01 


-2.18982846-04 
-2.1 799782E-04 
-2.16132756-04 
-2.l322928 c -04 
-2.0915779F-04 
-2.0 3 793 9 8 C — 04 
-1.97166966-04 
-1.89 330046-04 
-1.8042 3426-04 
-| . 70722 1 1 F— 04 
- 1 .6055 1 81 c -04 
-i. 50332986-04 
-1.40530 74^-04 
- 1 . 31 7733 5 c — 04 
-1.24619806-04 
- 1 • 196381 4 6— 04 
-1.1 7481 OC 6-04 
-1.1864379^-04 
-1.23543826-04 

- 1.32 548 19*=- 04 
-1.45808956-04 

- 1.6 1322 196— C4 
-1.84873346-04 
-2.10017^0^04 
-2.38187976-04 
-2.6845065F-04 
-2.998078?F-04 
-3.31161696-04 
-3.61037186-04 
-3.8815705 6-04 
-4.1121 82 5 c — 04 
-4.28753386-04 
-4.39656636-04 
-4.42949146-04 
-4.3 780606 E-C4 
-4.2381446 6-C* 
—4.00840 3 7 C -04 
-3.691 33 7 8 r -04 
-3. 292465 8 c -04 
-2.8225033F-04 
-2.29294 056—04 
- 1.72 104 8CF-04 
-1 • 12278 19 F— 04 
-5.18592 95F-05 

7.16677956-06 
6.3il3363 c -C5 
1. 1400342 E— 04 
1.58455266-04 
l .9511661 E— 04 

2. 232032 3E-04 
2.4219708F-06 
2.52081 506-0* 
2.53186786-04 
2.46352326-04 
2.32573116-04 
2.13348466-04 
1.90139426-04 
1.64 736386-04 
1.38688026-04 
1.13732226-04 
9. 11 709326-05 
7.22786536-05 



5.14025508-01 
5 .14338496-01 
5.1465140C-O1 
*>. 14964466—01 
5.152 77456-01 
5.15590448-0! 
5* 15903436—01 
5.16216426-01 
5.16*29418-01 
5.16842408-01 
5.17! 55398—01 
5.17468388-01 
5.1 7/81 3 78-01 
5.18094*58-01 
5.18407348-01 
5.1 8720338-01 
5. 1903332 c -0 1 
5.19366318-Oi 
5.19650305-01 
5 • 19S 72298-0 1 
5.20285285-01 
5.205982/ c -0l 
5.20911268-01 
5.2i224?5 c -01 
5.21537238-01 
5.21850228-01 
5.22163215-01 
5.22476205-01 
5.22 789195-01 
5.23102188-01 
5 • 234 | 51 7 C -01 
5.23728168-01 
5.24041158-01 
5.24354148-01 
5.24667128-31 
5.24980118-01 
5.25293108-01 
5.2560609^-01 
5.25919088-01 
5.262 12075-01 
5.26545068-01 
5.26858058-01 
5.271 7104E-01 
5.27484038-01 
5.2779702^-01 
5.28110005-01 
5.28422998-01 
5.28735988-01 
5.29048978-01 
5.29361965-01 
5 • 79674958—0 1 
5.29987945-01 
5.30 300935-01 
5.30613928-01 
5.30926918-01 
5. 31239896-01 
5.31552886-01 
5.31865878-01 
5.32178868-01 
5.32491 85E-01 
5*32804848-01 
5.33117838-01 


5 . 79071 40 E-05 
4.85954716-05 
4*45 766 1 0 5-06 
4.57248108-05 
5.15687996-05 
6.14776256-05 
7.4453564E-05 
8.95440548-05 
1.05542676-04 
1 .2 14561 OP -04 
1.36179176-04 
1. 4887087^-04 
1.58850116-04 
1.6558181 6—04 
1.66818046-04 
1.6851451 6-04 
1. 64803546-04 
i. 580 84448-04 
1.4881 383^—04 
1.37529618-04 
1.24923816-04 
1.1156056E-04 
9. 8 102 74fc £—05 
8.49922 3 O c -05 
7. 268I581E-05 
6.14477116-05 
5.14390186-05 
4*269161 4 £—05 
3.51 07575E-05 
2.04744156-05 
2.252 7896 c -05 
1 .69718316-05 
1. I4| 83?9 c -05 
5.55411626-06 
-9.69258528-07 
-8.39486378-06 
-1.69021866-05 
-2 .6633665^-05 
-3.75976216-05 
-4.97589576-05 
-6.29786938-05 
-7.70573566-05 
-9.1 7593078-05 
-1.06750176-04 
-1.2176302 c -04 
-1.36452? 3^-04 
- 1*505121 t c - 04 
-1.63679848-04 
-1.75707618-04 
-1.86423316-04 
-1.55675468-04 
-2 .0340 7C7 c -04 
-2.09oi237F-04 
-2. 14336476-04 
-2.17682206-04 
-2.19796026-04 
-2. 2085932 r - 04 
-2.21072766-04 
-2.20641386-04 
-2.19803688-04 
-2.187494 78-04 
-2.17669188-04 
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5.3343082F-01 
5.33 74381F-01 
5. 34 050806-01 
5 • 34 369796—0 1 
5. 34682 776—01 
5.3*995766-01 
5.3 C 08 75F-01 
^.3 >o 21746-01 
5 . 3593473 c -01 
5. *6247726-01 
5. 365607 1«=-01 
5.36073 706—01 
5.37186696-01 
5 • 3749968F-0 l 
5. 37812666-01 
5. 38 12565^-01 
5.38438646-01 
5 • 38 7 51 6 3* -01 
5 . 390 6 4 626-01 
5.3937761^—01 
5.39690606-01 
5 .4000*596— 01 
5.4O31658 r -01 
5 .406*9576-01 
5 • 4094256 c -0l 
5 .412^5546-01 
5.41 5 68 5 3 F— 01 
5.41 86152 E-0 l 
5.42194516-01 
5 .42 5 375 0 C -01 
5.42 820496-01 
5. 431334 06-01 
5 . 4*446476-01 
5.4375946^-01 
5.4^072456-01 
5. 44 3 854 36 -0 l 
5.4469842 c -01 
5.45011416-01 
5.4532440F-01 
5.456*7396-01 
^.45950386-01 
5.462633 76-01 
5.46 5763 6F -01 
5.46889356-01 
5.47202346-01 
5 • 475l533 c -0l 
5. 4782831 F-01 
5.401 41 30F— 01 
5.484 54296-01 
5 .48 767286— 01 
5 .49C802 7 c -0 1 
5 .49393266—01 
5 • 497062 5--0 1 
5 .50019246—01 
5.50332236-01 
5 • 506 45226-0 1 
5.50958206-01 
5. 512 71 1 9 C -01 
5.51584186-01 
5.518971 76-01 
5 • 522 10166-01 
5.525231 56—01 


-2.1671 '<*446—04 
-2.15976706-04 
-2.15541276-04 
-2.1540604^-04 
-2.15550536-04 
-2.15893766-04 
-2. 16296036-04 
-2.16611226-04 
-2.16622436-04 
-2.1614711 ^-04 
-2.1 49805 5 6— 04 
-2 • 129361 4F-C4 
-2. 0981759^-04 
-2.05520706-04 
-1.99991506-04 
- 1 .9325221 c - 04 
-1.85388346-04 
-1.76590016-04 
-1.67133016-04 
-1.57358 5 6 C — 04 
-1 .477707 76-04 
-1.38864066-04 
-1.3116956^-04 
-1.25359176-04 
-1.2 19933 7 C — 04 
-1.21593126-04 
- 1. 2473241 6-04 
-1.317545*6-04 
-1.4291633P-C4 
-1.583U58P-C4 
-1.7783803E-04 
-2.01207556-04 
-2.27884316-04 
-2.571078 76— C4 
-2. 88003036-04 
-3.19338596-04 
-3.49914636-04 
-3.70466136-04 
-4.03432456-04 
-4.235381 36- C4 
-4.37562066-04 
-4.44277346-04 
-4.42883306-04 
-4.32 704 1 5 r — C4 
-4.1348^116-04 
-3.85278926-04 
-3.48567036-04 
-3.04079455-04 
-2.53096436-04 
-1.97023006-04 
-1.37483436-04 
-7.65539276-05 
-1.60068135-05 
4. 19909456-05 
9.5761 50 5 e - C5 
1.4347350E-C4 
1.83858646-04 
2.1578626 c -04 
2.38644476-04 
2.52214166-04 
2.56712146-04 
2.52705726-04 



5. 5283614E— 01 
5.53149135-01 
5.53462125-01 
5. 53775116-01 
5.5408810 e — 01 
5.54401 08E-01 
5.54714075-01 
5. 5502 7065-01 
5.55340055-01 
5.55fc5304*-01 
5. 55966035—01 
5.56279026-01 
‘'.565920 IF— 01 
5. 56905 OOF -01 
5.5721 7996-01 
5.57530976-01 
5.57843965-01 
5. 581 5695F-01 
5 . 5846994F-0 1 
5. 5878293E— 01 
5. 59095925-01 
5 . 594C891 C -01 
5.5972190F-01 
5. 600348 9* -01 
5.t034788F-01 
5.6066087*— 01 
5.6097385^-01 
5.61286845-01 
5.fcl59583"-01 
5. 6191 282 E -01 
5.62225815-01 
5.62538805-01 
5.62851795-01 
5.631 64 78F-01 
5.63477775-01 
5.6379076**01 
5.64103745-01 
5.64416736-01 
5.64 725726-01 
5.65042716-01 
5.65355706-01 
5. 6566 86 9E— 01 
5.65981686-01 
5.6629467F-01 
5.66607666-01 
5.66920656-01 
5.67233646-01 
5.6 754662 c -01 
5.678596i: ; -01 
5.681 7?60 c -01 
5.68485596-01 
5.68798586-01 
5.65! 11575-01 
5 . 6942456*— 01 
5.69 73 755 c -Ol 
5. 70050545-01 
5.7036353**01 
5.70676516-01 
5. 7098950E-01 
5. 71302495-01 
5.71615485-01 
5. 719284 7*-0l 


2.412806 76-04 
2.23614686-04 
2.01407556-04 
1.76176235-04 
1.49803786-04 
1.23812876-04 
9.98718835-05 
7.42194326-05 
6.28449585-05 
5. 150528 7E-05 
4.54995346-05 
4.48308355-05 
4.91857946-05 
5.78604335-05 
7.0C760116-05 
8.47180535-05 
» .00769936-04 
1.17000525-04 
1.32430015-04 
1.46025655-04 
1.57074335-04 
1. 64994795-04 
1.69412995-04 
1.70244725-04 
1.67606925-04 
1.61756755-04 
1.53203’ 75-04 
1.42453265-04 
1.3020482 5-04 
1.17062035-04 
1.0360772*-04 
9.04419945-05 
7.79759075-05 
6.6 5669C05-05 
5.63684445-05 
4.74480576-05 
3.97302906-05 
3. 3C562C06-05 
2.71342686-05 
2.16452566-05 
1.62602765-05 
1.05786775-05 
4.3135023E-06 
-2.85932925-06 
-1. I14O495*-05 
-2.06485895-05 
-3. 14716315-05 
-4.35581825-05 
-5.68007355-05 
- 7. 1006041*-05 
-8.59125705-05 
-l.0I24272*-04 
-1.16667775-04 
-1.31871246-04 
- 1.465365 CE-04 
-1.60340475-04 
-1.73040925-04 
-1.84436135-04 
-1.94357025-04 
-2.02711495-04 
-2.09472325-04 
-2.146 70195-04 
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5 * 72241466-01 
5. 72 55445 c — 01 
5. 72867446-01 
5.73180436-01 
5 . 73493426-01 
5.73806406-01 
5 . 741 1939 c -0 1 
5.7443?38F-0i 
5 • 74745376-0 l 
5.75058366-01 
5. 7537135P-01 
5. 75684346—01 
5.75997336-01 
5. 7631 032 c -01 
5.76623316-01 
. 7693630^-01 
.7724928E-01 
.77562276-01 
. 776 75266-01 
. 7818825^-01 
. 785 01 24 c -0 1 
. 78814236-01 
.791 27?2=-0 1 
.79440216-01 
.79753206-01 
.00066196-01 
.803 79176-01 
• 80692 1 6 c -0 1 
.81005 ! 56-01 
. 0131 Bl4 r -0 1 
.81631136-01 
.81944126-01 
.82257116-01 
5.82 5701 0=-01 
5 • 82 8 3309=-01 
5.031 9608E-01 
5.8350907^-01 
5.83822056-01 
5.841 35046-01 
5.84448036-01 
5. 8476102 =-01 
5 • 050 7401 c -01 
5.85367006-01 
5 • 8569999 c -01 
5. 0601 2906-01 
5 • 86 3 2 597E— 0 1 
5 • 86638966 -01 
5 • 86 c 5 19 4 c -0 1 
5.9726493--01 
5.87577926-01 
5.8789091 = -0 1 
5.68203906-01 
5.88516896-01 
5.88629886-01 
5.89142 B 7E-01 
5.8945586F-01 
5. P976885F-01 
5 • 9000184^-01 
5.9039482^-01 
5 ♦ 907078 16-0 1 
5 . 91 020006-01 
5.91333796-01 


— 2.1838124E— 04 
-2.20747156-04 
-2.21927156-04 
-2.22Ill62 c -04 
-2.2 1551906-04 
-2.20458756-04 
-2.19054526-04 
-2.17525046-04 
-2. 160692 7 c -04 
-2.14814416-04 
-2. 1386466=— 04 
-2.13267146-04 
-2.13018086-04 
-2.13078706-04 
-2.13366016-04 
-2.1370792 c -04 
-2.13946526-04 
-2.1387764P-04 
-2.13334 006-04 
-2. 12062436-04 

- 2 . 0989300 = -04 
-2.06675636-04 
-2.0231 5lQ=-04 
-1. 96736636-04 
- I . 900 1 449=— 04 
-l. 82257226-04 
-1. 73655546-04 
-1.645475? c -04 
- 1 .553121 6 p -04 

- 1.4641530=-04 
-1.38362566-04 
-1.31 79998=— 04 
- I .27288966-04 
- 1.2536958=- 04 
-1.26659586-04 
-1. *157491^-04 
- 1 .40440 1 5 c -04 
-1.63489926-04 
-l. 7070955=- 04 

- 1 * 9 L 9 40 ? 1 6 - 04 
-2. 1674965=- 04 
-2.445311 3 C - 04 
-2.74518226-04 

- 3 .0557432=— 04 
-3. 36637356-04 
-3.6625025=-Q4 

- 3 .93081 3 56-C4 
-4. 1 582703=— 04 
-4.32991886-04 
-4.434331 76-04 
-4.4615376=- 04 
-4. 40293466-04 

- 4. 2543576 = - 04 
-4.01445406-04 
-3.68500416-04 
-3.2 7422536-04 
-2 * 79066 C6 6 -04 
-2 .247095 9= -04 
- 1.6615568=— 04 
-1.05043526-04 
-4.35839856-05 

1.64320496-05 



5.91646785-01 
5.91959778-01 
5.92272 768-01 
5.92585758-01 
5.92898 768-01 

5 *532 1 1 73F-01 
5.93524718-01 
5.93837 705-0 1 
5.94150698-01 
5.94463688-01 
5.94776678-01 
5.95089668-01 
5 .95*026 58—01 
5.95715648-01 
5.96028635-01 
5.96341628-01 
5.96654618-01 
5.96967598-01 
5.97280588-01 
5.9 75935 75-01 
5. 97906568-01 
5 .982 1955 c -01 
5.99532548-01 
5.98845538-01 
5.99158528-01 
5.59471518-01 
5.99784508-01 
6. 30097488-01 
6 . 004 1 0475—0 l 
6.0072^468-01 
o. 01036458-01 
6. 01349448-01 
6.01662438-01 

6 * Cl 975428-01 
6 .02288418-01 
6 • 02 60140 C — 01 
6. C291439E— 01 
6.03227388-01 
6.03540368-01 
6.03853358-01 
6 .0416634 c -01 
6.04479338-01 
6. 04792328-01 
6.05105318-01 
6.05418308-01 
6 *C5 731 298-01 
6.06044288-01 
6 .063572 78-01 
6.06670258-01 
6 .06983248— 0 1 
6.07296238-01 
6.0760922 c -01 
6 « 079222 1 8—0 1 
6.08235208-01 
6 * 085481 9 c -0 i 
6.08861188-01 
6.091 741 78-01 
6.C9487168-01 
6.09000158-01 
6 • 10 1 1 31 36-0L 
6.10426128-01 
6.10739118-01 


7.30031978-05 
1.24200268-04 
1.68606268-04 
2.04860528-04 
2.32190158-04 
2.50119845-04 
2.58657428-04 
2.582886 75-04 
2.49791978-04 
2.34470568-04 
2. 13742958-04 
1.89193075-04 
1.627281 0E-O4 
1.35938475-04 
1 • 106055 85-04 
8.8037063F-05 
6.95 I0657 c - C 5 
5*58001 70E-05 
4.74754688-05 
4.46127408-05 
4.69850798-05 
5 .4 1 20942 c - 05 
6.51457045-05 
7.91576098-05 
9.49461918-05 
1.11451365-04 
1 .2746821 E-04 
1.42019115-04 
1.54265865-04 
1.63500155-04 
1 .6932603 -—04 
i. 71530528-04 
1. 70140918-04 
1.65428695-04 
1.57766675-04 
1.47749118-04 
1.35941498-04 
1.23052698-04 
1.09712125-04 
9.64677195-05 
8.38442035-05 
7.21734828-05 
6.17140118-05 
5.25464878-05 
4.46485705-05 
3.78355875-05 
3.18982215-05 
2. 64 736095-05 
2.12296875-05 
1 • 58071 33 c -05 
9.81655815-06 
2.99901208-06 
-4.93885775-06 
- 1.4 1406 3 18-05 
-2.46834655-05 
-3.65786705-05 
- 4.9 7066 3 8 r - 05 
-6.39116098-05 
-7.89380338-05 
-9.45150995-05 
-1.10330638-04 
-1.26018275-04 



6.1105210E-01 
6.1136509F-01 
6.U67808E-01 
6.1199107E-01 
t. 12334 06E-01 
fc * 1261 705F-01 
6* 1293004F— 01 
6 • 13243026-01 
6. 1355601 F— 01 
6 • 13869006—01 
6 • 141 81 99F— 01 
6.1449498F-01 
6.148C797E-01 
6.15120966-01 
6* 15433956—01 
6* 1574694E— Ot 
6.16059936-01 
6.16372926-01 
6 • lt»68590 c -0i 
6# 1699889F— 01 
6.17315886-01 
6.1 7624876-01 
6. 1793786^-01 
6. 18250856-01 
6 • 18563 84F-01 
6.1b °^P-01 

6 • 191 89 1 01 

<.19502«l M 
< . 19815 796 'i 
6.2012878^-01 
6.204M 7 7<*-01 
6.2075476F-01 
6.2106775F-01 
6 • 21 380 746—01 
6.216937 3F— 0 1 
6.22006 72E-01 
< . 223 19 71 E-01 
6.22632 70 c -0l 
6 .22 94 5 6 9F— 01 
6 • 23 2586 7 C -01 
6.23571666-01 
6 • 23 88465F-01 
6 • 241 9 7 64F— 0 l 
6. 24510636-01 
6. 2482362 c -01 
6.251 366 1F-0 1 
6. 254496 OF-Ot 
6. 2576259F- 01 
6 • 2607 558F-0 1 
6 .263 8856E-01 
6 • 26701 55E-01 
6*2 7014546-01 
6.2732753E-01 
6 . 27640526-01 
6 • 27953 51E-01 
o • 2826650E-01 
6.285 7949E— 0 1 
6.28892486-01 
6 .292054 76-01 
6.2951846F-01 
6.2983144F-01 
6.3014443F-01 


-1.41259826-04 

- 1 • 5576 1246—04 
- 1*69191 836-04 
-1.81 34646F— 04 
-1*92 048136-04 
-2*011 54326—04 
-2*0 86208 1 c — 04 
-2.14437686-04 
-2. 186739 7F-04 
-2 • 2 143 669 E— 04 
-2*228731 3 p -04 

- 2 *232042 3F— 04 
-2*22626236-04 
-2.2138336F-04 
-2.1969391F-04 
-2. 1778204*^-04 
-2.158380 6E— 04 
-2.14053426-04 
-2. 125606 3F-04 
-2.11430306-04 
-2. 1070 76 3 C — 04 
-2. 10386546-04 
-2. 103966 3F— 04 
-2.10646586-04 
-2.10972486-04 
-2. 1 120254^-04 
-2.111483IF-04 
-2. 1057334^-04 
-2.09286946-04 
-2.07114536-04 

- -2.0391 86 5 r - 04 

-1.99581656-04 
-1.9409271^-04 
-1.8 7499196-04 

- I • 799386 3F— 04 
-1.71663746-04 
- 1 • 62 986 00 04 
-1.542 8046 6-04 
-1. 460959S*=-04 
-1.389534 7*^-04 
-1.33394766-04 
-1.30079736-04 
-1.29559186-04 
-1.323141 3F- 04 
-1.38 7652 76-04 
-1.49252016-04 
-1.6388116^-04 
-1.82617036-04 
-2.05178436-04 
-2.3113664E-04 
-2.59726036-04 
-2.900631 3 c -04 
-3.211236 76-04 
-3.5153172 c - 04 
-3.79998566-04 
-4.0517747^-04 
-4.25528366-04 
-4.3987322^-04 
-4.4693599F-04 
-4.45830C2F-04 
-4.15907346-04 
-4.16768946-04 
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6 .304 5742 6-0 i 
6.30770416-01 
6 • 31 08 340 F— 0 1 
6 . 31 39639 c — 01 
6.31 7003 86-01 
6.32022376-01 
6.32335366—01 
6 • 32648356-01 
6.32961336-01 
6.3327432E-01 
6 • 335873 1F-0 1 
f .33900306-01 
6. 342 1 3296-01 
6.34526286-01 
6 . 34 9 392 7 £—01 
6 • 351 5226F-01 
b. 35465256-01 
6.3577824^-01 
6.36091236-01 
6 • 364042 1 c -0 1 
6.36717206-01 
6.37030196-01 
6.37343186-01 
6.37656176-01 
6.37969l6 c -0l 
6.38? 82156-01 
6. 3859514^—01 
6. ”*89081 36—01 
6 .39221 12F— 01 
* .39534106-01 
6.3984709P-01 
6.4016008^-01 
6.40473076-01 
6. 4( ^8606^-01 
6.41099056-01 
6. 41412046-01 
6.41 725036-01 
6.42038026-Oi 
6.4235101^-01 
6.42663996-01 
6.42976986-01 
o. 43289976-01 
6.43602966-01 
6.43915956-01 
6.44‘>?F94r-0! 
6.4454193^-01 
6. 44854926-01 
6.4516791 r -0 1 
6.45480906-01 
6.4579?89E-01 
6.461 C6 8 7E-01 
6.4641^866-01 
6.46 73285P-01 
6.4 704 5 84f -01 
6.4735883^-01 
6.476 71826-01 
6.47984816-01 
t . 4829 780E-0 1 
6.48610796-01 
6.4892378C-01 
6.49236766-01 
6 • 49549 75f-0 l 


-3.08507856-04 

-3.51526796-04 

- 3. 0662 293 F— 04 
-2.548983 7P-04 
-1.9796850E-04 
-i.3742582 c -04 

- 7. 5397695 E-05 
-1.37171696-05 

4. 530046 2 E— 05 
1 .00146706—04 
1.48700006-04 
1.8966243E-04 
2.2 1981 52F-04 
2. 449 20*06—04 
2.5828605E-04 
2.6 231836E-04 
2.57621266-04 
2.45334906-04 
2. 26 71447 E— 04 
2.C352057E-04 
1. 7734288E-04 
1.50130046-04 
l .23500566-04 
9.91307326-05 
7.82489496-05 
6.1972416F-05 
5.0916739P-05 
4.54232576-05 
4. 54330 l 86—05 
5.0 50645 8E- 05 
5.99619986-05 
7. 2911694 c -05 
8.81 91 3076-05 
I.0472295F-04 
1.21238136-04 
1.36741316-04 
1.50199126-04 
1.60915846-04 
1.66306 1 5 C — 04 
1.72097156-04 
1.72235236-04 
1.68859356-04 
1.62369316-04 
1.53213116-04 
1.42083276-04 
1.29594836-04 
1.16377976-04 
1 • C 3 1 1 400 p - 04 
9.02 775 7 7F — G c 
7. 03 399 8 5 C — C c 
6. 75*06606-05 
5 • 8G4l 3 33 c -0* 
4.98405 3 7 C -C C 
4.2*51754F-C5 
3 . 680526 5 C - 06 
3. 14240546—0*' 
2 .63 1 4965F- C6 
2.11248956-0* 
1.55048776-05 
9.0690019^-06 
1.58535476-06 
-7. 177II89P-06 



6 .4986274F— 01 
6.5017573F-01 
6. 50488726-01 
6.5080171F-01 
6*511 1470^-01 
6.5142769E-01 
6.51740686-01 
6*52053676-01 
6*52366666-01 
6 • 526 7964 c — 01 
6.5299263E-01 
6 • 53305626—01 
6* 5361 06 IF- 01 
6*53931606-01 
6.5424459E-01 
6. 545 57586-01 
6.5487C57 r -01 
6* 551 8356F-01 
6* 55696 55 F— 01 
6*55809536-01 
6 • 561 22 52F-01 
6. 5643551^-01 
6.56748506-01 
6.5706149F-01 
6*57? 74486—01 
6.57687476-01 
6 • 580QC466-01 
6. 59313456-31 
6.58626449-01 
6.58939436-01 
6. 592 5241 F-01 
6* 595 5540F-01 
6. 59678396-0 1 
6*6019138^-01 
6. 60504376-01 
6.6081 736 c -0l 
6.61 l 3035F-01 
'• 6144334 6-01 
6-617 56336-01 
6.62069326-01 
6.62302306-01 
6.62695296-01 
6. 630082 8 F— 01 
6.633212 7F-01 
6.636 34266-01 
6.6354725E-01 
6.64260246-01 
6.64573236-01 
6.64886226-01 
6.65199216-01 
6.65512206-01 
6.65925! 06-01 
6.661381 7 r -01 
6.66451166-01 
6.6676415^-0! 
6.67077146-01 
6.6739013 F-0 1 
6.67703126-01 
6.68016116-01 


-1*7325391 F— 05 
-2*88664736—05 
-4*17388106-05 
-5.58063216-05 

- 7*085 1 7 10 c -05 
-8.65774876-05 
-1.02678906-04 
-1.18822146-04 
-1.34627036-04 
-1 .49805 1 7 c -04 
-1.64020776-04 
-1. 77021036-04 
-1. 8858402E-04 

- 1 • 98 563 86 p - 04 
-2.06872336-04 
-2.13471306-04 
-2.18402556-04 
-2*21 742256-04 
-2*2 36200 8 04 
-2.24244626-04 
-2.2 3808 1 1 6—04 
-2.2253818 — 04 
-2.2067598^-04 
-2.1845576 c -04 
-2.16097736-04 
-2.1 380076^-04 
- 2* 1 1 731 60 E- 04 
-2.1001660^04 
-2.08 734 50 c - 04 
-2.0793120 c -04 
-2.0754997^-04 

- 2 .0754540 E- 04 
-2.07791346-04 
-2.08149976-04 
-2.08454096-04 
-2.08444926-04 
-2.07929846-04 
-2.0672709C-04 
-2.04657626-04 
-2.01565656-04 

- 1.9 13 72 9 3 C -04 
- I .9 206 76 5 E— 04 
-1.657066^6-04 
-1. 78453996-04 
- 1 . 70 5 5860 F— 04 
-1.6232592 6-04 
-1.54234356-04 
-1 .4675525 c - 04 
-1.40396346-04 
- 1* 3581 771 c -04 
-1.33582086-04 
-l. 34193^8^-04 

- I • 3821522F-04 
-1.4600610^-04 
-1.57818886-04 
- 1. 7374 L 49 6-04 
- 1. 9366090 c -04 
-2-1 72941 1 F— 04 
-2.44022626-04 
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Input for Sample Case 2 


SAMPLE CASE 2 - VH « 0. 

40ATA ALPHA = .114* RADIUS * 4.3* NBLAOE * 4, PHI * 26., RPM * 35S.» 
VH = D.Ot X * 30. ♦ 50. » 10. • * = 15* BETA * .25* IPLOT * 1* 

COEM-1 » 1.. .46238, .05370, .0501 2, .03273* .01905, .00691* .01585. 
.02661* .00708, .00794, .00631, .01776, .00531, .00501, .00467, 34*0., 
C"'tt"P2 ~ 0., 40*, -25. , 20., 15., —30., 44*0.4 


Output for Sample Case 2 


ft OAT A 

ALPHA * 0.1L4E +00, 


ppTf r o.rse+oo. 


C1EFFI 

= 3.1E*01, 0 

.46238E+00, 

0.537E-01, 

0. 5012F— 01 , 

0.3273E-01, 


0.19055—01* 

0.891E-02, 

0.1585E-01 

, 0.266IP-0I 

, 0. 7C8F-C2, 


0.794P-02, 

0.63l*-02, 

O.l 778E-01 , 

0.5? IE— 02, 

0.501E-02, 


0.46 7E-0?. 

O.Ot 0.0, 

0*0, 0*0 • 

0.0, 0.0, 

0.0, 0.0, 0.0. 


0.0, 0«0« 

0*0, 0*0, 

0*0, 0.0, 

0.0, 0.0, 

C.O, O.C, 0.0' 


0»0* 0 *0 • 

0*0f 0 . 0 * 

0*0* 0.0, 
0.0, 

0*0, o.c. 

0.0, 0.0, 

0.0, O.C, 0.0 

CCEFF2 

= 0.0, 0.4E+02, -0.25E+0 2, 0.2E+02, 

0.15E+02, - 

0.3E+02, 0.0, 


0*0* 0.0, 

0*0, 0.0, 

0.0, 0.0, 

0. C , 0.0, 

0.0, 0.0, 0.0 


0.0, 0 .0 « 

0*0, 0.0, 

0.0, 0.0, 

0.0, 0.0, 

C.O, 0.0, 0.0 

IPLOT 

N 

OaCi o.o, 

0* Ot 0 *0 , 

* l, 

* 15, 

0.0, 0*0, 

0* 0 • 0.0, 

0*0, 0.0, 

0.0, 0.0, 

0. 0 , 0.0, 

0.0, 0.0, 

0.0, O.C, 0.0. 
C.O, 


NBl AOE - 4, 


PHI * 0.26F+02, 

RADIUS * 0.43F+01, 

PPM * 0. 355E+C3 « 

VH * 0*0* 

X = 0* 3E +02 * 0.5E+02, 0*1 F02t 

SEND 
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1 


1 


1 


SA*»Pl* CASE 2 - VH * 0. 


TIME 

INTERPOLATED PRESSURE 

l .67546626-01 

-9.8168699E-05 

1 .67859616-01 

-e.9946558E-05 

1.681 7260^-01 

-8.47941 74«=-05 

1 .68485596-01 

-8.28038576-05 

1.68 790586-01 

-8.4442538=-05 

1.6911157F-01 

-8.9974628E-05 

1.6942456F-01 

- 9.85602326-05 

1 . 69 73 755 r -0i 

-1.13229466-04 

I • 70050546—01 

-1.308201 IE-04 

l.7036353*-0l 

-1.5200C4OC-O4 

1 . 7067651 6-01 

-1. 76292546-04 

I • 70989506—01 

-2.0309C27E— 04 

1.71 30249E-01 

-2.31487006-04 

1.71*15406-01 

-2.6060^746-04 

1.7192 84 7 c -0 1 

-2.894 71 8 3 c — 04 

1.72241466-01 

-3.16870796-04 

1.72554456-01 

-3.417278I e -04 

1. 72 86 7446-01 

-3.629S606E-C4 

1.7318043^-01 

-3.7958169^-04 

1. I349342 c -01 

-3.9C57819E— 04 

I • 73806406-01 

- 3 .85243 44 F— 04 

1.741 1939 c -0l 

- 3. 9 30245 3 P - 04 

1.744*2306-01 

— 3.835487C6— C4 

1 • 74745376—01 

-3.6675843 c -04 

1. 75058366-01 

-3.42 7802 IE-04 

1.75371356-01 

— 3.12131666— C4 

1.7568434 r -01 

-2.7S48I03E-C4 

I .75°97336-01 

-2.33675? 8F-C4 

1 • 7631 032 p -01 

- 1.8798304^-04 

1. 76623316-01 

-l. 3963477E— 04 

1.76936306-01 

-8.99262076-05 

1.7724928^-01 

-4.04621 786-C5 

1.77562276-01 

7.39359236-06 

1.77875266-01 

5.2359301E-05 

1 . 7818825F-01 

9.30738386-05 

1 .78501246-01 

1.2853553^—04 

1 • 788 14236—01 

1.5 795 2726-04 

1.79127226-01 

1.80703P 7E— 04 

l • 794402 1 6—01 

1.56530376-04 

1 . 79753206-01 

2 .05404106— 04 

1.80066196-01 

2. C 76449 36- 04 

1.80379176-01 

2.03 77306E— 04 

1.80692166-01 

1.94502026— 04 

l.SlOOSlSF-O’ 

1.8C9C751E-04 

1.81310146-01 

1.64020566-04 

1.01631136-01 

1.4 49462 3E- 04 

1.81944126-01 

1.25006516-04 

1 .82257116-01 

1.05300956-04 
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I 


) 

l 

i 

» i 


1.8257010**01 
1.82883096-01 
1.8319608 c *0l 
l .835C9076-01 
1 . 83822056— 01 
l • 84! 35046-01 
1.84448036-01 
1.8476l02 c -0l 
1.85074016*01 
1.8538700^*01 
1.85699996-01 
1.8601298^01 
i • 863259 7 C *01 
1 . 86 638966*01 
l • 86 95 l 94- *01 
1.87264936-01 
I .8757792**01 
1 .67890916*01 
1 . 88203906—0 1 

l .8351639^^01 

i .88825886*01 
1 - 89142876*0 1 
l .89*5586 -—01 
I . 89 7688 5 r * 01 
1.9008164^-01 
I .5039482F-01 
i -99 7078 1 c *0l 
1.51 020806— 0 1 
1.9131379^*01 
1-91^4678**01 
1. 91955776-01 
1.922 7276*^-01 
1.92585756*01 
1.92698746*01 
1-93211 73**01 
1-5352471 c -01 
1 -938 37705*01 
1.941 50696*01 
1.94463686*01 
1 -94 7766 76*01 
1 .95089666*01 
1.9540265—01 
I .95715646*01 
l • 96026o3 c *0l 
1 .96341626-01 
1 -9665461 F*01 
1 .56967596-01 
1.97280586-01 
1 -975935 7 C *01 
1.97906566*01 
1 • 982 I955 P *0I 
1*985? 254F-01 
i • 98845536—01 
1 .99 158526-01 
1.99471518-01 
1 .99 78450^-01 
2.00097486-01 


8.68646586*05 
7.01982266— 05 
5.73809406-05 

4- 76296586— C5 
4. 1684601 c *C5 
3.96315426-05 
4.13138296-05 
4. 6396805^-05 

5- 44101486*05 

6- 4641043^*05 
7.63895196-05 
8.89224106-05 
i. 01 358076-04 
1.13015256*04 
1.23256156-04 
1.31468246-04 
1.37278996*04 
1.403943 86*04 
1.4C 700296*04 
1. 382 091 9 C - 04 
1.33102 866*04 
1.25638746-04 
1.16171986-04 
1. 05162 C96-04 
9.304O664F-Q5 
8.02670646-05 
6.73124256-05 
5.45615126-05 
4. 23136846*05 
3.0855206 c -05 
2. 03542966-05 
1.08905586-05 
2.46503866-06 

-5.C0757I 66*06 
-1.16 536366*05 
-1.76 36003 F-05 
-2 .32014 74 c *05 

- 2 *85502536-05 
*3. 3889169 c *05 
-3.94537866-05 
*4.53923906-05 
*5.18294886*05 
*5.88437776-05 
*6.6474597^*05 
*7.46881506*05 
*8.34436426-05 
-9.26285606-05 
- 1*02 1145 3 C — 04 
*1.11 742 96 6— 04 
-1.21374486-04 

- I • 3082049 c *04 
-1.39962696-04 
-1.48671976-04 

1.5685S92E-C4 
-1. 6442681 F-04 
-1. 7135340^-04 
-1.77637156-04 


TO 



< ! 1 1 1 

I 

• ♦ 

; ' i ' 

$ i 

* f » i i i 

! 

t 


2. 00410476-01 

-1.83290936-04 

2 .00723466-01 

-1.86419746-04 

2.01036456-01 

-1.93045536-04 

2.01 34944 c -01 

-1.97268856-04 

2.01662436-91 

-2.01177176-04 

2.01 97542 e -0l 

-2.04856895-04 

2.02208416-01 

-2.03365136-04 

2.02601406-01 

-2. 11750536-04 

2.0291439E-01 

-2.15044586-04 

2.03227386-01 

-2. 19219236-04 

2.03540366-01 

-2.2 1227036-04 

2.03853356— 01 

-2.2398429E-04 

2.04166345-01 

-2.2 6358436-04 

2.04479336-01 

-2.28219686-04 

2.04792325-01 

-2.29410305-04 

2.05105316-01 

-2.29707366-04 

2.0 5418306-01 

-2.23950616-04 

2.05731296-01 

-2.26952466-04 

2. 060442 8F-01 

-2.2355 7696-04 

2.063572 7*-01 

-2.13686766-04 

2.C66 7025E— 01 

-2.12241316-04 

2.0698324F-01 

-2.04228806-04 

2.07296236-01 

-1.94719126-04 

7.07609276-01 

-1.8 3847786-04 

2.07922216-01 

-1.71024546-04 

2. 08235206-01 

-1.54887986-04 

2 .085481 9 c -0 1 

-1.45700636-04 

2.0886118F-01 

-1.32427926-04 

2.09174176-01 

-1.18664436-04 

2.09487166-01 

-l.C 792 3 3 76-04 

2.09800156-01 

-9. 78442 746-05 

2.10113136-01 

-8.89630986-05 

2.10426126-01 

-8.4 7941 756- 05 

2.10739116-01 

-8-28038576-05 

7. 11052106-01 

-8.44425896-05 

2.1I365C96-01 

-8.99746276-05 

2. 1167808 F-01 

-9.9560230E-05 

2.11991076-01 

-1.13229486-04 

2.12304066-01 

-1.30820116-04 

2.1261705E-01 

-1.52000406-04 

2.12930046-01 

-i.7629254=-04 

2.13243026-01 

-2.03090266-04 

2.13556016-01 

-2.31486996-04 

2.13869006-01 

-2.60603746-04 

2. 14181996-01 

-2.89471826-04 




Error Diagnostics 

The following table gives the error diagnostics as the message and the 
reason for getting the message. All of the errors are fatal. 


Message 

TIME HAS NOT REACHED TMIN FOR BLADE I 


MACH NUMBER IS NOT LESS THAN ONE. 
PROGRAM ABORTED. 

BLADE I FAILED TO REACH TMAX IN ^00 
ITERATIONS FOR INTERPOLATION SET K. 


ITH ELEMENT OUT OF ORDER. NO INTER- 
POLATION. 

IORDER = 0 

ONE ELEMENT IN X ARRAY 
INSUFFICIENT NUMBER OF POINTS. NO 
INTERPOLATION. 


Reason 

The ith blade has been rotated back 
1000 time increments, but its initial 
arrival time has not yet reached a 
close enough range to TMIN. 

This analysis is to be applied only 
to velocities less than Mach 1 * 

The ith blade has net reached TMAX 
after 200 points have been computed 
starting from its respective initial 
time; occurred during interpolation 
set k. 

> Error messages from IUNI 
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APPENDIX 
SUBROUTINE IUNI 


LANGUAGE: FORTRAN 

PURPOSE: To interpolate a univariate function using conventional first 

or second order Lagrangian interpolation. The routine will 
perform multiple table look-up for a set of functions defined 
over the same set of independent parameter points ♦ The 
points in the independent parameter array car be unequally 
space 1, 

USE: CALL IUNI (NMAX,N\X,TTAB,Y ,IORDER,XO ,Y0 ,IFT,IFRR) 

NXAX An input integer specifying the maximum 

first dimension of the functional value 
array as given in the dimension state- 
ment of the calling program. 

?l An ^nput integer specifying the actual 

number of independent parameter points, 
where N £ NMAX. 

X An input array dimensioned at least N 

in the calling program. This array con- 
tains the values of the independent 
parameter. The elements of the X array 
mi be strictly monotonic. 

N7AB An i .put integer specifying the number 

of dt pendent variable tables. 

Y A two-dimensional input array whose 

columns contain the dependent variable 
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tables. The Y array should have first 
dimension NMAX and second dimension at 
least NTAB. 

An input integer specifying the order of 
interpolation to be used 

- 0 Zero order interpolation: The 

first value in the dependent vari- 
able table is assigned as the inter- 
polated val:.e of the function. 

= 1 First order interpolation. 

* 2 Second order interpolation. 

The point at which interpolation is to 
be performed. 

An output array containing the inter- 
polated value of each function. This 
array should be dimensioned NTAB in the 
calling program. 

An input/output integer parameter with 
the following functions 

INPUT: Initialize routine IUNI nd 

check monotonicity of the X array. 

- -i Whenever a new X array is input, 

this value of IPT must be specified 
by the user so that a monotonicity 
check of the X array will be 
performed. 



iuni 


OUTPUT: INDEX PCI?rTEF 

* k Indicates that x. < x <x . On 

k- 0 x+1 

the next call, the previous value 
of IPT is used as a pointer to 
begin the search for the subinterval 
containing the next interpolation 
point . 

Whenever the point x 0 is not contained 
in the interval delimited by the X 
array, extrapolation is performed to 
estimate the function value. In this 
case the value of IPT is returned as: 

= 0 Indicates x^ < x'l) if the X 

array is in increasing order, or 
Xq > x(l) if the X array is arranged 
in decreasing order. 

= N Indicates > x'N) if the X array 
is in increasing order, or 
Xq < x(N) is the X array is arranged 
in decreasing order. 

IERR Integer error parameter generated by the 

routine. 

= 0 Normal return 

= J The J-th element of the X array was 
out of order. No interpolation 
performed. 
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IUNI 


RESTRICTIONS:' 


METHOD: 


= -1 Zero order interpolation performed 
because I0RD2R = 0. 

= -2 Zero order interpolation performed 
because only one point was in the 
X array. 

= -3 Insufficient number of points 

supplied for second order inter- 
polation. No interpolation 
performed . 

= -U Extrapolation was performed. 

The parameter I ERF should be tested 

return to the calling program. 

The points in the X array must be arranged in strictly increasing 
or decreasing order of magritude. 

Whenever a new X array is provided, the parameter IPT - -1 must 
be input by the user to insure the routine will be properly 
initialized and the monotonicity of the X array checked. 

In table look-up for scientific computing, consecutive inter- 
polation points tend to fall in the same region of the indepen- 
dent variable table. The parameter IPT provides an index to 
begin the search for the subinterval containing XC. IPT is 
used to initialize the routine on the first call; thereafter 
it preserves the location in the table of the previous answer. 
This feature was implemented to minimize the amount of search 
time. 
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IUHI 


ACCURACY: 

REFERENCES: 

STORAGE: 

SUBPROGRAMS 

USED: 

FORTRAN 

FUNCTIONS: 


If Y, , Y_, . . . ,Y„ are the values of the function evaluated 
12 S 

at the tabulated values x 1 * x 2 f **** ] Si of the independent 
variable, then the interpolation equations are as follows: 


Second Order: The index i is chosen so that 

ho - x i I * k - i i«l 18 8 


^ X i+2* X i^ 


- V : W*i - x 0 ))( V2 - x o ) 


(*i+l ’ *i } 




- V - 


w»i+i ~ x o ))(x i : x p ) 


u 


i-2 


“ X i+l J 


First Order: For <, x Q <, x^ +1 


, . ( v i ~ V - ( .y. v 

0 1 \*i ~ v 


Accuracy is a function of the order of interpolation used. 


Nielson, K. L., Methods in Numerical Analysis; pp. 87-91. 
Ralston, A., A First Course in Simerical Analv is ; pp. U2-U6 

U04 octal locations 


None 


SIGN 
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SUBROUTINE PSEUDO 


LANGUAGE: 

PURPOSE: 

USE: 


COMPASS 

To create and write an appropriately named Plot Vector File. 
Through linkages set up by an initial call to PSEUDO, all 
subsequent graphics data generated by the user will be routed 
through one of the PSEUDO entry points and written on the Plot 
Vector File. 

CALL PSEUDO 
or 

CALL PSEUDO (FN) 

FN file name left- justified with zero fill. Default 

file name is SAVPLT. 


Example : 

CALL PSEUDO 

This will establish a Plot Vector File named SAVPLT. 
CALL PSEUDO (6IMYFILE) 

This will establish a Plot Vector File named MYFIIE. 

NOTE : The Plot Vector File (or Files) will usually be 
written to disk (as opposed to tape) and may be 
postprocessed following user program termination via 
appropriate specification of one or more PLOT control 



PSEUDO 

RESTRICTIONS: 


METHOD: 


STORAGE: 

SUBPROGRAMS 

USED: 


(1) An initializing call to PSEUDO tvith or without a 
file name argument) must be made prior to any calls 
to CALPLT or any other graphics output routine. 

(2) Every Plot Vector File should be terminated with a 999 
pen code, CALL CALPLT (0.0,0.0,999)* The transmission 
of the 999 code will cause an EOF write on the Plot 
Vector File, and the file will temporarily be closed. 
Thus, any given Plot Vector File will contain only 
one 999 pen code and/or one EOF. 

(3) To continue plotting execution following + - ansmission 
of a 999 code to a current Plot Vector File, the user 
program must call the PSEUDO processor to create new 
Plot Vector File (i.e., CALL PSEUDO ( 6LMYFIL2 ) ) . 

In addition to entry PSEUDO, this processor contains two other 
entry points, namely PLT9999 and PLT999® • An initializing call 
to PSEUDO will set PLT9999 into the processor switching mech- 
anism (PLOTSW) . Subsequent plot dat \ generation will then be 
routed via CALPLT, PLOTSW, and PLT9999 and written on the Plot 
Vector File. The entry PLT9990 is used record special 
purpose data from routines NFRAME and PLTSTOP. 

2155 g locations total for direct subprograms 

NUMARG, PLOTSW 
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SUBROUTINE CALPLT 

LANGUAGE: FORTRAN 

PURPOSE: To move the plotter pen to a new location with pen up or down. 

USE: CALL CALPLT (X,Y,IPEK) 

where 

X,Y are the floating point values for pen 

movement . 

IPEN = 2 pen down 

= 3 pen up 

Negative IPEN will assign X a 0, Y * 0 
as the location of the pen after moving 
the X,Y (create a new reference point). 

= 999 Writes a terminating block address of 

999 to terminate the Plot Vector File 
^ and all further processing is skipped. 

CALL CALPLT (0.0,0. 0,999) 

RESTRICTIONS: All X and Y coordinates must be expressed as floating pcirt 

inches (actual p-tge dimensions) in deflection from the origin. 

A 1J2L ABATING BLOCK APDRE0S (IPEN 999) M**ST BE GIVEN AS THE 
LAST PLOTTING INSTRUCTION 1EFCRE ENDING A PROGRAM WHICH OSES 
A1.Y OF THE PLOTTER SUHR , J , J rir ”S- THIS IS TO BE SURE THAT ALL 
PLOTTER 1I.oTRUCTIC 73 APE WTiTTEN ON THE PLOTTER TAPE. 
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CALPLT 

METHOD: 


STORAGE: 

SUBPROGRAMS 

USED: 


The main subroutine in the graphics language is the CALPLT 
subroutine. All other special purpose subroutines eventually 
call CALPLT either directly or indirectly. This routine moves 
the pen in a straight line between the present pen position 
and another pen location to which the prograaer wishes the pen 
to be moved. 

In order to cause such instructions to be written, the programer 
specifies the coordinates of the point to which the pen is to 
be moved and whether the pen is to he moved in a raised or 
lowered position. This is accomplished by the FQRTRAH 
instruction: 


CALL CALPLT (X,Y, IPEH) 

A new plot origin within a frame is created by a call to CALPLT 
with the argument, IPEH, negative. The CALPLT routine then 
moves the pen to X,Y; stores this location as (0,0), that is, 
a new origin. Each call to CALPLT with IPEH negative will 
establish a new plot origin within frame. After a frame of 
plotting is completed , a call to HFRAME should be used to 
establish a new frame origin. This new frame origin is 
referred from the last frame origin, not from the last plot 
origin created by a call to CALPLT with the argument , IPEH, 
negative. 

251g locations 6000 series 


PLQTSW, STRCALL, LOCATE 
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LANGUAGE: 

PURPOSE: 

USE: 


SUBROUTINE AXES 


FORTRAN 

To draw a line, annotate the value of the variable at specified 
intervals with or without tic narks , and provide an 
identification lalei. 


CALL AXES (X,Y, THETA, S,0RG,SFX,TMAJ,TMIH,3CD, HEIGHT, NOCHAK) 


where 

X,Y are the coordinates in floating point 

inches of the starting point of the 
axis with reference to the plotting 
area origin as established by CALPLT. 

THETA is the angle of rotation measured 

counter-clockwise from the X-axis in 
floating point degrees. Normally, 

THETA is 0° for an X-axis and 90° for a 
Y-axis. 

S is the length of the axis in floating 

point inches. Should be a multiple of 
THAJ. 

♦S will generate tic marks. 

-S will eliminate tic marks. 

ORG is the functional value to be assigned 

to the origin the value of the 

first scale) in floating point. 
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SFX is the adjusted scale factor for the 

array to be plotted (change in value 
per inch). 

MOTE: Values of ORG and SFX which will 

produce a reasonable scale may be 
calculated usi v* ’♦'he subroutine ASCALE 
or BSCALE. 

TMAJ is the distance in floating point inches 

for major tic marks (0.25 inches high). 
Numbers are placed on the axis the 
major tic marks in accordance *ith the 
values of ORG and SFX. The numbers 
written along the axis are adjusted 
♦: be between lOOC.OO and 0.01 in magni- 
tude. Imediately after the last 
number on the axis is placed the caption 
xlO eX ^ 9 where exp is the required 
exponent . 

If the values are integer multiples, the 
decimal point and decimal places are 
eliminated. A negative TMAJ will cause 
the actual value to be written instead 
of the adjusted value. 

UCN is the number of divisions per inch in 

floating point for minor tic marks 
(0.125 inches high). To eliminate 
minor tic marks the following may be 
used: 


TMIN = 0. 
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RESTRICTIONS : 

STORAGE: 

SUBPROGRAMS 

USED: 


BCD is the character label for the axis 

(see NOTATE routine). 

HEIGHT is the height of the full-size 

characters in tne BCD title. Numbers 
at the tic marks will be (O.T5*HEIGHT) 
high. HEIGHT is in floating point 
inches. 

If HEIGHT * 0., all annotation will be 
eliminated. 

NOCHAR is an integer specifying the number of 

chaiacters in BCD title. A negative 
NOCHAR places the annotation on the 
clockwise side cf the axis and a 
positive NOCKaR places the annotation 
on the counterclockwise side of the axis. 
NOCHAR * 0 is net allowed. If it is 
desired to have nc label, then the BCD 
parameter should be 1H , and NOCJAP 

eitber *1 or -1. 

Only perpendicular axes are recoomendec . 

1^51g locetions 6000 series 

CALPLT, NOTATE, NUMBER, ROUND, SIN, COS, V^ERE 
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LANGUAGE: 

PURPOSE: 

USE: 


SUBROUTINE LINE 


FORTRAN 


To drew a continuous line through and/or draw a symbol at each 
successive data point (stored in an array). 

CALL LINE(XARRAY*IARRAY,N,K t J,L s HEIGHT) 


where 


XARPAY, 

YARRAY 


N 


K 


J 


are the names of arrays containing the X 
values and Y values , respectively, to 
be plotted. Values must be in floating 
point. 

is the number cf points to be plotted. 

is the interleave factor which speci- 
fies the sequence in which data ore 
stored. 

s 1 indicates that values are stored 
sequentially. 

* 2 indicates that values are stored 
in every other location in the 
array, etc. 

is positive for line and symbol plot, 
negat ive for symbol only j lot . The 
Mgr.itude specifies the alternate 
briber of data points at which to 
aiot a symbol. 
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RESTRICTIONS: 

STORAGE: 

SUBPROGRAMS 
USED: - 


* 0 for line plot. 

* 1 for symbol for every data point. 

= 2 for symbol for every other data 
point , etc. 

L is an integer describing symbol to be 

used; see NOTATE routine for list. 

HEIGHT is the desired symbol height in float- 

ing point (see ROTATE routine). 

LINE expects the adjusted minimuas and scale factors as 
described in ASCALE. 

352g locations 6000 series 

CALF1T, ROTATE, WHERE 
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SUBROUTINE NFRAME 


LANGUAGE: FORTRAN 

PURPOSE: To provide users specific means of executing frame advance 

movements on any plotter device via an appropriate frame 

oriented device postprocessor. Frame advance distances are 
generally defined to be incremental from current frame origin 
(i.e., comparable to frame advance executions for the DDI or 
252 CRT devices). CALL NFRAME is intended to be used as a 
frame advance mechanism, not as a plot origin offset. 

USE: CALL NFRAME 

or 

CALL NFRAKE(H,V) 


where 

H and V are Horizontal 'parallel to device X) 

and Vertical 'parallel to device Y ) 
distances from the current frame origin. 
H and/or V must be expressed in float- 
ing point inches. 

The short form CALL NFRAME will cause the device postprocessor 
to execute a fra re advance move parallel to the device X 
{horizontal) axis. The movement will be the maximum horizontal 
distance used in inches plus h incnes, where h will be an 
increment appropriate tc the particular device {0 < h £ 2"). 

The frame advance will he rounded to whole ^nches. 

When H an'* V parameters are provided or. the NFRAM* call, 
only f Growing values are permissi v le: 
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CALL NFRAME(H,0. ) Frame advance K" horizontal 
CALL NFRAME(0. ,V) Frame advance V" vertical 
CALL NFRAME(H,V) Frame advance E" hy V" 

CALL NFRAME(0. ,0. ) Return to current frame origin 
CALL NFRAME'H,-V) Frame advance K" by -V" 

This should be used to execute a frame advance move parallel 
to horizontal axis and establish a new origin for roll paper 
plotters and return to origin for flatbed plotters. 


RESTRICTIONS: 


1. This routine is intended for use only in concert with the 
frame-dependent graphics postprocessors. 


2. This routine must be used in any case vhicb may require 
AUTO modification of Plot Vector File data by a graphics 
postprocessor. 

3. The frame advance distances specified by E and/or V 
should always be at least slightly greater them the 
intended usable frame size. H and V are ignored by the 
Varian and DDI postprocessors. 

U. The H dimensions of a frame advance may not be negative 
(H ^ 0). For purposes of frame stacking, V may be either 
negative, z®ro, or positive. 


STORAGE: ?6g memory word^, oJC 6000 series 

SUBPROGRAMS NUMARG, CALPLT, PlT99y8, ABORT 
USED: 
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